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PERIODIC ORBITS ABOUT AN OBLATE SPHEROID 

BY 

WILLIAM DUNCAN MACMILLAN 

§ 1. Introduction. 

The orbit of a particle about an oblate spheroid is not in general closed 
geometrically. The motion of the particle is not, therefore, in general, periodic 
from a geometric point of view. But if we consider the orbits as described by 
the particle in a revolving meridian plane which passes constantly through the 
particle several classes of closed orbits can be found in which the motion is 
periodic. The failure of these orbits to close in space arises from the incom- 
mensurability of the period of rotation of the line of nodes with the period of 
motion in the revolving plane. When these periods happen to be commensur- 
able the orbits are closed in space and the motion is therefore periodic, though 
the period may be very great. Indeed, it seems that most of the difficulty in 
giving mathematical expressions for the orbits about an oblate spheroid rests 
upon the question of incommensurability of periods. The difficulty arising from 
the motion of the node can be overcome by the use of the revolving plane, but 
other incommensurabilities, such as that introduced by the eccentricity, can not 
be eliminated in this manner. 

Orbits closed in the revolving plane are considered most conveniently in two 
general classes : (1) Those which reenter after one revolution, (2) those which 
reenter after many revolutions. The existence of both classes is established in 
this paper and convenient methods for constructing the solutions are given. 
Orbits which reenter after the first revolution are naturally the simpler and will 
be considered in the first part of the paper. Those lying in the equatorial 
plane of the spheroid become straight lines in the revolving plane, and it is 
shown that within the realm of convergence of the series employed all orbits in 
the equatorial plane are periodic. When the orbits do not lie in the equatorial 
plane there exists one, and only one, orbit for any arbitrarily assigned values of 
the inclination and the mean distance. These orbits reduce to circles with the 
vanishing of the oblateness of the spheroid. 

In considering orbits which reenter only after many revolutions the differ- 
ential equations are found to be very complex, and one would despair of proving 
the existence of these orbits by the ordinary direct computation of the necessary 
coefficients. However, a proof of their existence and a method for the construc- 
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56 W. D. MACMILLAN : PERIODIC ORBITS [January 

tion of the solutions are given by the aid of certain theorems, which are here 
established, on the character of the solutions of non-homogeneous linear differ- 
ential equations with periodic coefficients. 

These periodic orbits of many revolutions involve five constants, four of which 
are entirely arbitrary, and the fifth subject only to the limitation that it shall 
satisfy certain commensurability conditions. One constant, only, is missing for a 
complete integration of the differential equations. These orbits are all symmetric 
with respect to the equatorial plane. 

§ 2. The differential equations. 
The differential equations of motion of a particle about an oblate spheroid are * 

x- + y 2 — \z 2 



d 2 x 

dY 2 = 




d 2 y 
dt 2 ~ 


it 3 [ + » 


drz 
df = 


R i A+ 10° 



l_sv 
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The symbols employed are defined as follows : x, y, z are rectangular coordi- 
nates, the origin being at the center of the spheroid and the a;y-plane the plane 
of the equator, k is the Gaussian constant, If is the mass of the spheroid, b is 
the polar radius of the spheroid, ft, is the eccentricity of the spheroid, 



, , -o k 2 M\^ b 2 x 2 + y 2 - 
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we obtain one integral of areas, namely, 




dy dx 

^ x dt- y dt=^ 





■]• 



That is, the projection of the area described by the radius vector upon the equa- 
torial plane is proportional to the time. We have also the vis viva integral 

For further integration we are compelled to resort to series. Poincare has 
shown in his Les methodes nouvelles de la m.ecanique celeste that if certain 

* Moulton, Celestial Mechanics, p. 1 13. 
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conditions are fulfilled it is possible to obtain periodic orbits represented by 
power series in a parameter when periodic orbits are known for zero values of 
this parameter. These conditions are at least partially fulfilled in the present 
problem, for the right members of the differential equations are analytic in 
x, y , z and the parameter fi, . Furthermore, for /u. = , the equations reduce 
to the ordinary two-body problem for which periodic solutions are known. It 
is our purpose to show that the remaining necessary conditions also are fulfilled 
and that periodic solutions persist for values of /x =)= . These periodic series 
are very satisfactory, for the general character of the orbits represented by them 
is easily obtained. The solutions thus found are rigorous, but they are not 
general, their existence depending upon special initial conditions. 

It will be advantageous to transform the differential equations to cylindrical 
coordinates by the substitutions 

x = ar cos v, k 2 3f=?i 2 a i , 

y = ar sin v , c. = ck V Ma , 

(4) 

z = aq , tit = t , 

R = aV^+q\ l 6a l=0l 
After these substitutions equations (1) become 

(a) /'-™ 2 = p-- ?2? [l+p^y^>=+...], 

(5) (b) rv" + 2rV = 0, 

the accents denoting the derivatives with respect to t. 

The integral of (5) is r^v ' = c, by means of which v can be eliminated from 
the first of these equations. After the elimination the equations take the form 



r 3 — 4rq 2 



( c ) W ' = J- 

The first two of these equations are independent of the third so that r and q 
may be considered as rectangular coordinates in a revolving plane which passes 
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always through the z-axis and through the particle itself. The problem is thus 
reduced to the consideration of the motion in this plane, for, when r is known, 
v is obtained by a simple quadrature from (6c). 

§ 3. Surfaces of zero velocity. 
The velocity integral in the revolving plane is 

, 2 , 2 _ 2 2 r 2 -2q 2 2 c 2 

r + 1 = (rTjr^y* + 3 (V +Y) 1 ' ^ + ^ 2 + ° 2 ' 

If we put the velocity equal to zero the resulting equation represents a two- 
parameter family of curves. For assigned values of the parameters c and c 2 
there is defined a certain curve in the revolving plane. On one side of this 
curve the motion is real while on the other side it is imaginary. For values of 
c 2 < this curve is closed, and the motion is real on the inside. As the plane 
revolves this curve generates a surface of the general form of an anchor ring. 
For /a 2 = this curve belongs to the ordinary two-body problem, and its equa- 
tion is 

2 

W+<ff 



vi- l_ ^2 \i — Jl + C 2 _ " ' 



The motion is elliptic, parabolic, or hyperbolic according as c 2 is negative, zero, 
or positive. Putting 

r = p cos cf> , q = p sin <f> , 

we find, on solving for p. 



^=^[- 1± \( 1 + cos^]- 



For negative values of c 2 this equation represents two closed ovals which do not 
enclose the origin. If c 2 c 2 = — 1 the ovals shrink upon the points p = — 1 /c 2 , 
<j> = and 7r . The corresponding orbit is therefore a circle in the equatorial 
plane. As c 2 approaches zero the ovals open out rapidly and approach the 
limiting curves 

c 2 _ 
P== 2~eoa r 4>' 

For values of c 2 > there is but one positive value for p, which is 

If c 2 4= none of these curves crosses the axis <f> = 90°. But if c = we have 
the circle p = — 2/c 2 inside of which the motion is real when c 2 is negative. 
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For values of p? 4= , but sufficiently small, we can put 

r = (p + p)cos<f>, q = (p + p)sm<f>, 

and solve for p as a power series in p? . We find 

1 2-3cos 2 d> , 

which is the correction to be applied to the corresponding curves in the two- 
body problem. 

Part I. 

ORBITS WHICH REENTER AFTER ONE REVOLUTION. 

§ 4. Symmetry. 

Returning to the differential equations (6 a) and (66) we observe that if we 
change 

rinto + »', q into — q, t into — t, 

the differential equations remain unchanged. Hence, if at some epoch, t = t , 

r = a, r = 0, 
<7=0, ? '=/3, 

that is, if at the epoch t = t , the particle crosses the r-axis perpendicularly, it 
follows from the form of the differential equations that the orbit is symmetrical 
with respect to the r-axis and with respect to the epoch t = t . In other words 
r will be an even series in (t — t ) and q will be an odd series in (t — t q ). If 
now at some other epoch, t =t 4- T, the particle again crosses the r-axis per- 
pendicularly the orbit will be symmetrical with respect to this epoch also. It is 
clear therefore that the orbit is a closed one, and that the motion in it is peri- 
odic, for, by hypothesis, at t = t — T it must have been at the same point and 
moving with the same velocity in the same direction. The motion is therefore 
periodic with the period 2 T . Hence with these initial values sufficient condi- 
tions for periodicity are that at t = t + T 

r = q = . 

From the integral of area, t/= c/r 2 , it follows that if r is periodic v will have the 
form 

v = constant x t + periodic terms. 



c 2 
(«) r " = r s- 


1 


0> 2 ^M 4 


(6) «'=£• 
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§ 5. Existence of periodic orbits in the equatorial plane. 
If q = equations (6) reduce to 



(8) 



The first of these equations is independent of the second and can be inte- 
grated separately. It represents motion in a straight line in the revolving 
plane. It admits the constant solution 

r = i, ci = i + e\^ + e\^ + --., 

which represents a point in the revolving plane, or a circle in the equatorial 
plane. In order to investigate the oscillations about this point let us put 

r = 1 -f ep , c ! =Cj+ ec, 

where p is a variable whose initial value is arbitrarily assigned, e is an arbitrary 
parameter corresponding to the eccentricity in the two-body problem, and e is a 
parameter to be determined so that p shall be periodic. 

On substituting these values (8a) and expanding as power series in e, the 
terms independent of e cancel out and it is possible to divide through by e. 
The equation then becomes 

p " + [1 _ 0> 2 - 30> 4 + • • •] /> = e [1 - 3/se + 6/aV - IO/jV + • •] 

+ p 2 e[B - 40J/* a - 150> 4 + ...] 

(9) + />V[ - 6 + 100* p* 4- 466>~ /*' + • . • ] 

+ /3 4 e 3 [10 - 20<?> 2 - 1110> 4 + • • • ] 



We can simplify this equation somewhat by dividing through by the coefficient 
of p in the left member and then substituting 



The equation then becomes 

(10) 5 + / 3 = 8 C 1 - 3 / 3e + 6 ' 32e2 - 10 P 3 e' , + ---] + [34-a 1 ]p 2 e 

+ [ - 6 + o 2 ] P 3 e 2 + [ 10 + aj />V + 
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where 

a 2 = + 40> 2 + (40* + 2801 ) /*' + ■ ■ ■ » 

« 3 = - 1O0> 2 — (lOtfJ + 810 2 ) ^ + • • • . 

Equation (10) can be integrated for p as a power series in 8 and e . Let us 
take the initial values 

P=-l, p' = 0. 

By Poincare's extension of Cauchy's theorem we know that the solution of equa- 
tion (10) having the prescribed initial values exists and converges provided 8 
and e are sufficiently small, for all values of T in the interval = TS T, where 
T is finite, but otherwise arbitrary. The condition for periodicity is simply 

(11) p'=0atT=r. 

If we choose T= ir an inspection of equation (10) shows that for e = the 
solution for p is periodic with the period 2ir whatever may be the value of 8, so 
that equation (11) must carry e as a factor. After integrating equation (10) 
we find that the condition (11) is 

(12) = — [| + a, J7rSe — [ § + \a x + -*. 2 a\ + |« 2 ] ire 2 + higher degree terms. 

Aside from the factor e there remains an equation in which the linear terms in 
e and 8 are present, and which may be solved. We find 

(13) S = (_l + ...) e +.... 

If this value of 8 be substituted in equation (10) it will then admit periodic 
solutions for p as power series in e having the period 2tt for all values of e suf- 
ficiently small. Furthermore the solution as a power series in e with the pre- 
scribed initial conditions is unique. 

§ 6. Existence of periodic orbits not lying in the equatorial plane. 
For /u. 2 = the differentia] equations (6) admit the circular solution 

r = 1, </ = 0, v=r, 

(14) r f c 2 =l. 
/ = , q = , v = 1 , 

In order to investigate the existence of orbits not lying in the equatorial plane 
but having the period 2tt for /Lt 2 =(= , let us put 

(15) r=l + p, ,7 = + 0-, c 2 =l + e, 
and take the initial conditions 

p = a, p = 0, o- = 0, o-' = /3//,. 
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The conditions for periodicity (§ 4) are then 

/a' = <7 = at T=7T. 

We have three arbitrary constants at our disposal, a, /3, and e, and two condi- 
tions to be satisfied. Hence one constant remains undetermined. We will 
therefore let /3 remain arbitrary and determine a and e so as to satisfy the 
two conditions. After making the substitutions (15) and expanding, equations 
(6) become 

p" + P = e - Bpe + 3p 2 + f <r 2 - 8\ p 2 + Qp 2 e - 6p 3 - Qp<r 2 

(16) 4 \pd\p? 4 higher degree terms, 
a-" 4 a = 2>pcr — 6/a V 4 f /> 3 — 0"# 2 /* 2 + higher degree terms. 

In order to integrate these equations let us put 

(17) (i+j + k>0), 

The /o r; . and a\. 4 can be found by successive integrations, the constants of 
integration being determined so as to satisfy the initial conditions. In the series 
thus obtained take t = ir . The two conditions for periodicity give us the two 
equations : 

(a) p! r =„ = = a l €?+a 2 ea+a s <?+a i (?a+a 5 ea 2 +a e a 3 +a 7 ep 2 +a 8 afjt. 2 +a 9 iM l -\ , 

(6) <r T= „ = = /S/i [V + ^e 2 + 6 3 « 2 + 6 4 e« + 6 5 ^ 2 + . • •], 

where the a ; and 6 ; are constants which have been found from the actual 
integrations to be distinct from zero. Equation (18a) involves only the even 
powers of p while (186) involves only the odd powers. After dividing (186) 
by /S/x, we can solve it for e as a power series in a and p 2 of the form 

(19) e = c x a 2 + c 2 p? + c 3 a 3 + c t ap? + c 5 /x 4 + • • • . 
On substituting (19) in (18a) we obtain a series of the form 

(20) (a) = d 1 ap, 2 + d 2 a 3 + d 3 p* + d 4 a 2 it 2 + d 5 a* + • ■ ■ . 
If in this equation we make the substitution 

we obtain 

(c) = /**[/ l7 +/,/** 4 /.7M 1 +•••], 
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which can be solved uniquely for 7 as a power series in p? . This solution sub- 
stituted in (206) gives a as a power series in /a 2 , and this value of a substituted 
in (19) gives e as a power series in p? . We thus have a solution 

a = p 2 P l (p?), e = /a 2 P 2 (/a 2 ), /3 arbitrary. 

Newton's parallelogram shows that equation (20«) has two additional solutions, 
but as they are imaginary we do not stop to develop them. 

§7. Existence of periodic orbits in a meridian plane. 

If in equations (6) we take the area constant c 2 = the motion of the particle 
is in a meridian plane ; that is, the revolving plane has ceased to revolve and the 
orbit in this plane is the true orbit. After changing to polar coordinates by 
the substitution 

r = p cos <f> , o=p sin <f> , c 2 = , 



the differential equations become 

(21) 



« ^ , 1 - f + f cos 2<ft + * cos 4^ , 



» . A sin 2<f> — 4 sin 4<f> ,„ „ 

For p? = we have the periodic solution 

P = 1 • # = t , 

that is, a circle. For /a 2 =(= let us introduce p and o- by 

p = l -t p, (f> = t + a , 

with the initial values 

p = a , |o' = , o- = , o-' = /3, 

where a and /3 are two new arbitraries. By Poincare's theorem /), //, <r and 
cr' are expansible as power series in a , /3 and /t 2 with t entering the coefficients. 
The conditions for periodicity are that at t = it 

p' = a- = . 

If we perform the integration and then set t = 7r , we obtain the two following 
conditional equations : 



(a) o- T=lr = = a 1 a + a 2 /3 + a 3 a 2 + a 4 ot/3 + a 6 /3 2 + a 6 /t 2 + • • • , 

(b) ,w = 0= & 3 a 2 + & 4 «/3 + & 5 /3 2 +0V+& 7 /i'+..., 
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where the a. and b. are constants which have been found from the actual inte- 

t % 

grations to be distinct from zero. 

The first of these two equations can be solved for a as a power series in /3 and 
y? . This expression for a substituted in (b) gives rise to an equation of the 
form 

(c) = c,/V + c 2 /3' 3 + c 3 /3V + c if j. 4 + ■■■. 

This equation has the same form as (20a) and can be solved in the same way, 
giving a unique real value for /3 as a power series in p? vanishing with p? . 
This expression for /3 substituted in the equation for a gives a as a power series 
in /j. 2 , vanishing with p? . Therefore real periodic orbits exist for /la 2 =|= 
which are analytic continuations of circular orbits for p? = . 

We have thus proved the existence of the three following classes of periodic 
orbits which have the period 2tt , the generating orbits being circles : 
I. Orbits lying in the equatorial plane ; 
II. Orbits not lying in the equatorial plane; 

III. Orbits in a meridian plane. 

§ 8. Construction of periodic solutions in tlie equatorial plane. 

Let us first consider orbits in the equatorial plane. We retake the differ- 
ential equations (8) and by means of the transformations there given we pass at 
once to equation (10), which is explicitly 

fp + P =8[l-Bpe+6p 2 e 2 + • • •] + [3-30>»-(0f +60» )/*<+■ • > 2 e 

( 23 ) + [-6 + W\p? + (4&i + 286;) p' + ■ • ->V 

+ [1O-1O0> 2 -(1O^ + 810 2 )^+...>V+.... 

It was shown in equation (13) that 8 can be expanded uniquely as a power 
series in e in such a manner that the solution for p as a power series in e will be 
periodic with the period 2ir . Since the series is periodic with the same period 
for all values of e sufficiently small, it follows that the coefficient of each power 
of e is itself periodic. Since the solution exists and is unique, it must be pos- 
sible to determine the 8 uniquely by the condition that the solution is periodic. 
In the existence proof it was shown that 8 vanishes with e . Therefore p and S 
have the form 

(24) p ==Po + Pi e + p 2 e 2 + p 3 e* + ..., 8 = 8^ + 8.^+8^ + -... 

The p. are to be determined by the integration of equation (23) and by the 
initial values 

dp 
P = -l, ^ = 0atT = 0. 

The 8 are to be determined in such a manner that the p. shall be periodic. 
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Substituting (24) in (23) and equating the coefficients, we find : 
(«) d r+Po=0, 

(c) d *£ + P . i = 8 2 -Zp l) 8 l + [3-Z0lS-(0* + 66l)S + ...]2p l)Pl 
(25) ai 

+ [_6+40> 2 + (40 4 + 280 2 > 4 +...>3, 



(<*) fj' + p„ = 8-3 / > o S„_ 1 +[3-30> 2 -(0*+60 2 ) / , 4 + •••]2/ V >„_, 



These equations can be integrated in succession. The solution of (a) which 
satisfies the initial conditions is 

(26) p = — cos T . 

Since the initial conditions are independent of e, every p. except p must vanish 
at T = 0. Substituting (26) in (256) and integrating, we have 

(27) ^=8^1-008 T) + [3-30> 2 -(0 4 +60|> 4 +. • •] [J- J cos T-i cos 2T] . 

The constants of integration in equation (27) have been determined so as to 
satisfy the initial conditions, but the constant 8^ is as yet undetermined. 
Substituting these values of p and p { in (25c), we find 

^ 2 2 + P 2 = [8 2 + 8,(3-30> 2 +•••) + 3 -60> 2 + ...] 

+ [8,(-3 + 60> 2 + { 20? + 120 2 } M 4 + •••) 

(28) + ( - 3 + 120>* - { -V-<?i - 96> 2 } /i 4 + • • •] cos T 
+ [8,(3 - 30> 2 + ...) + 3 - 180> 2 + •••] cos2T 

+ [3 -40 2 ,x 2 + ••■] cos 3T. 

In order that the solution of this equation may be periodic the coefficient of 
cosT must be zero. This is the condition which determines 8 i and consequently 

8 l = - 1 + 2&\^ + (|0f - 0f> 4 + . ■ -. 

Trans. Am. Math. Soc. 5 



66 W. D. MACMIL.LAN : PERIODIC OKBITS [January 

With this value of B y equation (28) becomes 

f^l + l° 2 = P,+ 3^ 2 +- • •] + [- M\f + . ■ •] cos 2T + [3-40^ 2 + ■ • •] cos 3T. 

The solution of this equation which satisfies the initial conditions is 

ft =8 t (l-cosT) + [381^+...] + [t-y4>* + ...]conT 

+ [30> 2 + ...]cos2T + [- | + !«*/*» "+...]oob ST. 

The constant 8 2 is as yet undetermined. It is determined by the periodicity 
condition for p 3 in the same manner that B l was determined by the periodicity 
condition of p 2 . Without giving the details of the computation its value is 
found to be 

8, --60? /** + .... 

This method of integration can be carried as far as may be desired. In 
orderto show this let us suppose that p , • • •, p n _ l have been computed and all 
the constants determined with the exception of S n _ 1 . From (25«!) we have 

(30) d ^ + P n = B n -Bp S n _ l + lS-Mi^ + ...]2p p n _ 1+ f n ( Po , ..-,/.._,). 

Here ,f„(p , ■ ■ ■, P„_ 2 ) * s a polynomial in the p. and contains only known terms. 
p n _ l depends upon S n _ l in the following way : 

p n _ 1 = S„_j(l — cosT) + known terms. 

Equation (30) can therefore be written 

^ + P. = 8.+ [3-8^*+...]«^ 1 +[-3 + 6(??^ + ...]S |l _ 1 oo8T 

+ [3 — 30? /t 2 -+■ ■ • •] cos 2T + known terms. 

In order that the solution of this equation shall be periodic the coefficient of 
cos T must be zero. This condition determines S nl . The equation can then 
be integrated and the constants of integration will be determined by the condi- 
tions that 

d Pn 
dT 



p n = -f> = at T = . 



Everything is then determined with the exception of S n , and we have then 
p n = ( 1 — cos T ) 8 n -\- known terms. 
Substituting the values of & Y and 8 2 in the solution as far as it has been com- 
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puted, we find 

P = — cosT, 

ft = [1 + J^M 2 + (0 4 - 4flJ)^- ••] + [- 0> 2 + (- f* 4 + 30*>< + • • ■] cos T 

+ [-4 + |0> 2 +tt^ + 2 > 4 -"]cos2T, 
/> 2 = [- 30> 2 + • • •] + [I - |0> 2 + • • -]cob T + [30> 2 + • • ]cos 2T 

+ [-[t + ^> 2 +---]cos3T, 
S 1 = -l + 20> 2 + (f0J-0 2 >4 + ---, 
S 2 =-60> 2 + .... 

From these expressions we have the following series for r = 1 + ep: 
{a) r = 1 - e cos T + { [i + J^M 2 + (0[ - 40 2 > 4 + • • • ] 
+ [_0 2 /i 2 + (_7^ + 30 2 >*+...]cosT 
(31) + [_i + i0^ + (i0* + 03) At 4 + ...]cos2T}e 2 

+ { [ - 30 2 M 2 + • • • ] + [ | - i#> 2 + • • • ] cos T 
+ [30> 2 + • • • ] cos 2T + [- f + |0> 2 + • • • ] cos 3T } e 3 + • ■ ■ . 
Substituting this value of r in the equation (86) 



-£ 



r 2 dv . _ X _ c S + ee 



and integrating, we find 

(b) v -v = {[i + eis + m + mi)^ + ■■■] 

+ [*>*+(- f*f + V-W + • • •] « 2 + • • • } T 

(31) + {[2 + 20> 2 + (0 4 + 2 > 4 + ---] sinT } e 

+ { [20> 2 + (J^J - 69DS + • ■ ■ ] sin T 

+ [| + !#> 2 + ( - M 4 + |^ 2 )^ + ■ • • ] sin 2T } e 2 + • • -. 

Equations (31a) and (316) are the periodic solutions sought. If we return to 
the symbols defined in the original differential equations by means of equations 
(4), with the additional notation 

«l/l-0> 2 -30 2 V + --- = v, 
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we have the following expressions for the polar coordinates at any time t : 

R = a { 1 — e cos vt + [ J — J cos 2vt ] e 2 + [ f cos vt — f cos Bvt ] e 3 + • • • 



(32) 



+ ^ ^ 2 [ ( 2 % - T 3 o cos vt + fo cos 2rt ) e 2 + ( - ^ - ^ C os vt 
+ -j 9 cos 2vt + j»j cos 3^)e 3 + • • • ] + -j /x 4 [...]+...} , 



i 2 
•u — v = vt + 2 sin itf • e + | sin 2itf • e 2 H + - ^[(^ + T 8 T e 2 + ■ • • )itf 

(B3) " ° &4 

+ (| sin i4)e 2 +(§ sin vt + ^ sin 2i4)e 2 + • • •] + ~ 4 M 4 [- ••]+•• ' 



« 



These equations contain four arbitrary constants, a, e, v and t u . Since the 
differential equations of motion in the equatorial plane were of the fourth order 
these series, within the realm of their convergence, represent the general solu- 
tion. The expression for the radius vector, It , is always periodic with the period 
2<tt\v. At the expiration of this period v has increased by the quantity 

^[^ "'(ft + lV 2 +•••) + ■••] = 2^ 

in excess of 2tt; that is, the line of apsides has rotated forwarded by this 
amount. If © is commensurable with unity the orbit is eventually closed geo- 
metrically. If © = 1 1 J, where / and J are integers relatively prime, then, at 
t = 2Jrrjv, v = 2(1 + J)ir , and the particle is at its initial position with its 
initial components of velocity. The particle has completed I+J revolutions 
and the line of apsides has completed /revolutions. The mean sidereal period is 

This formula for the rotation of the line of apsides has an interesting appli- 
cation in the case of Jupiter's fifth satellite. On the hypothesis that Jupiter is 
homogeneous and taking its equatorial diameter as 90,190 miles and its polar 
diameter as 84,570 miles, the mean distance of the satellite as 112,500 miles, 
the eccentricity of its orbit as .006 and its mean sidereal period as ll h m 23 s , 
the above formula gives for the rotation of the line of apsides 1440° per year. 
The values derived from observations are somewhat discordant but are in the 
neighborhood of 880° per year. If we still keep the hypothesis that Jupiter 
is homogeneous in density and of the same oblateness as before, we are com- 
pelled, in order to relieve the discrepancy between theory and observation, to 
suppose that the value adopted for its polar radius was about 9,000 miles too 
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great. From the large reduction required it is clear that the hypothesis of 
homogeneity is very much in error. 

§ 9. Construction of periodic solutions not lying in the equatorial plane. 

By means of the area integral the problem has been reduced to the two 

equations 

c 2 r r 3 — 4irq 2 2 

r = r 3 - ( ,.2 + ?2)3 - (y + q y *i * + • ■ ■ « 

( 35 ) 

1 (r 2 + ? 2 )i (r^+YJi ^ + "" 

After the solution of these equations has been obtained the third coordinate is 
found from the equation 

c 

We have already proved [equations (14) to (20)] the existence of periodic solu- 
tions of these equations of the following type : 

r^\ + p 2 ^ + Pi ix i + ..., 
(36) q = q l n + q g fi 3 + q^ 6 + • • • , 

C 2 =l + C 2/ Lt 2 + C 4 ^ + 

with the initial conditions 

r'(0) = y(0) = 0, q'(Q) = fo, 

/3 being arbitrary. We can therefore integrate equations (35) so as to satisfy 

these initial conditions and determine the c. in such a manner as to render the 

j 

solution periodic. 

Substituting (36) in (35) gives the equations 

= [P* +P-2- kl - C 2 + *?]/*' + IX + Pi~ 3/>2 - 3?,? 8 + 6p 2 q[ 

(87) + W + ( 3c * - w ?)ft - Wrf - "J/* 4 + • ■ •> 

= ^ i+ q ^ + [ q '*' +q *~ 3 ^i- i?i+3% 1 ]m 3 + [?r+?5- 3 M 3 - i ?& 

38) +6^-3^+^^+^1+3^.-16(9^-^??]^+ • • •• 

Equation (37) contains only even powers of \x while equation (38) contains only 
odd powers. For the integration we have : 
Coefficient of p . 

(39) ' ' ?;' + ?, = o. 
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The solution of this equation satisfying the initial conditions is 

q t = /3 sin t . 
Coefficient of /j? . 

rf + ft- Itf + *-*?, 

= (i /3 3 + c 2 - 6\) - 4 /3 2 cos 2t. 

The solution of this equation which satisfies the assigned initial conditions is 

p 2 = (f/3 2 + c 2 - 2 ) + a 2 cos t + i/3 2 cos 2t. 
The constant c 2 will be determined by the periodicity condition on q 3 to be 

c 2 = 20 2 -/3 2 , 
and a 2 will be determined by the periodicity condition on p t to be 

a 2 =0. 
If we anticipate these determinations we have 

/ > 2 = (0 2 -i/3 2 )+i/3 2 cos2T. 
Coefficient of /* 3 . 

(41) 

= ( 3/3 2 + 3c 2 - 60 2 ) /3 sin t + |a 2 /3 sin 2t . 

In order that the solution shall be periodic it is necessary that the coefficient of 
sin t shall be zero. Therefore 

c 2 = 20 2 -/3 2 . 

Substituting this value and integrating, we find 

q 3 = /3 3 sin t — j« 2 /8 sin 2t. 
From the initial conditions we must have q' 3 (0) = and therefore 

But it will be shown in the next step that a 2 = , and consequently 

Coefficient of /a 4 . 
(42) p'; + p 4 = 3/> 2 + 3 ?1?3 - 6^ - -V41 + (40? - 3c 2 )/> 2 + J£0J ? ; + c 4 . 

Before expanding the right member of this equation we will avoid useless labor 
by first examining the coefficient of cos t which we know must be zero from the 
periodicity condition. We notice first that terms in the coefficient of cos t can 
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arise only through terms involving p 2 and ^ g , and second that all such terms 
carry a 2 as a factor. Since no other arbitrary enters the coefficient, we must 
take a 2 = . It can be shown by induction that the constant of integration a 2i 
(the coefficient of cos t) which arises in the integration of p 2 . is determined by 
the periodicity condition on p 2(i+0 , and further that its value is always zero. 
The proof is omitted for the sake of brevity. 

Substituting the value a 2 = in p 2 and q z and expanding the right member 
of (42), we find 

Pi + P< = K + °\ + V- i & ~ &P ] + C l 6 l & ~ T«P ] eos 2t + 11& cos 4t • 
The solution of this equation satisfying the initial conditions and the properties 
just mentioned is 

p i= [c 4 + B\ + ^.6\&--i^-\ + [- r V??/3 2 + ^/8«]oos 2t- &/8*oos4t. 
Anticipating the value of c 4 found below, we have 

p 4 = [ — 3^^ - y 0?/3 2 - & /3 4 ] + [ - j\d*f3* + T y3 2 ] cos 2t - J f /3 4 cos 4r. 
Coefficient of p? . 

?." + ?5 = b ftj2 3 + f ?I?S - 6 />2?i + 3 ?4?l - ¥ft?I ~ W 

(43) -8(?J?, + 16(?J/. f?1 + ^5rJ, 

= [ 8e 4 + 120* + 116> 2 /3 2 ] /3 sin T - 2 /3 3 sin 3t . 
From the periodicity condition we have 

c 4 = -40*--y-0 2 /3 2 . 
Integrating and imposing the initial conditions, we find 
q r> = - f 0\ /3 3 sin t + \8\ /3 3 sin 3t . 

This is sufficient to make evident the general character of the series. The r- 
equation contains only even multiples of t , and the ^-equation contains only odd 
multiples. The r-equation contains only even powers of p. and of t. The q- 
equation is odd in both these respects. The series are therefore triply even and 
odd. 

Collecting the various coefficients we have the following series : 

(a) r=l+ [(0 2 -i/3 2 )+i/3 2 cos 2t> 2 + [{-^\-\\6\^-^) 

+ ( - tV e \ F + 16 &) cos 2t - e^ 4 cos 4t] ^ + ■ • • , 
(6) ? =[^sinT]^+[O]^+[-|0 2 /3 3 sinT+^ 2 iS 3 sin3T] M 5 +---, 

(c) «-« =[l_^y_(|^ + i^ / 3 2 )^+...]T+[-i/3 2 sin2T]^ 

+ [( \\6\ /3 2 - \&) sin 2t + J^/3 1 sin 4t] p* + • • ■ , 
(<Z) c 2 = l + [2^- y 8 2 ]^ + [-4^--V-^/8 2 ]M 4 +--- 
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This solution contains four arbitrary constants, a, @,v u and t„. As is shown by 
equation (44c) the nodes regress, the rate of regression being 

[0> J + (§*? + *<?? /8>« +••■]• 

The generating orbit of these solutions is a circle in the equatorial plane. A 
circle having any assigned inclination might have been used, e. g., 



(45) r == i/l — s 2 sin 2 t, q = ssiar, v = tan _1 [i/l — s 2 tan t], 

where s is the sine of the inclination. The solution thus obtained would have 
been identical with (44). If we expand (45) as power series in s and then put 
s = /3/u , it will be found that the terms thus obtained are identical with the terms 
independent of 9\ in the solution which has been worked out. It might there- 
fore be of assistance in the physical interpretation of the constants to put /3fi = s 
in the series (44). 

§ 10. Construction of periodic solutions in a meridian plane. 

When the constant c is zero the motion is in a meridian plane. The equa- 
tions of motion are (21) 

» ,-2,1 -f + fooB2$ + icos4A 

p -p<t> + p i = - - -— «- 0; /*- + • • • , 

(46) 

a" , 9 'a.' £ sin 2<£ - £ sin 4<£ 
p4> + 2p <f> = -j 0;> 2 + • • -. 

We have already shown the existence of periodic solutions of these equations as 
power series in /u. 2 , which for /j, = reduce to the circle p = 1 , 4> — t . Let us 
put then 

/> = 1 + P^+p^ + ■■-, <k = t + <£ 2 M 2 + 4>X + • ■ -. 

Substituting these expressions in (46), expanding and collecting the coefficients 
of the various powers of /x , we find : 

= [ p ;'_3p 2 -2^;-f 6\+ \6\ cos 2t+ |0 2 cos4r> 2 + [p';_ 3^-2^-2^; 

-0; 2 + 3^ 2 + (3-6 cos 2t-cos4t)6' 2 p 2 -(-(-3 sin 2T-sin4r)6';</> 2 ] / Li 4 -(-- • •, 
(47) 

0= O;'+2/>; + i0J sin 2t-|0 2 sin 4t] M 2 + [ft + 2/>; + # ft +2/>#; 

+ (— 2 sin 2t + sin 4t)0 2 p 2 + (cos 2t — cos 4t)0 2 <£ 2 ] /*' H . 

The initial conditions are 

p'(O) = 0(O) = O. 
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Proceeding to the integration we have : 
Coefficient of ft 2 . 

(a) p.; - z lh - 2& = \e\ - %e\ cos 2t - \e\ cos 4t , 

(6) <^;' + 2 i >;=-i^sin2T+i^sin4T. 

By integrating (5) once we have 

(c) <f> 2 = - 2p 2 + \Q\ cos 2t - -^e\ cos 4t + c, . 
Substituting this value of cj> 2 in (a), we find 

(d) p' 2 ' +p. 2 ={2c 1 + %6\ ) — 6\ cos 2t - §6\ cos 4t . 
The integration of this equation gives 

(e) p 2 = (2c, + \9\ ) + c 2 sin t + c, cos t + \0\ cos 2t + -fo6\ cos 4t. 

Since |)j = at t = we must take c 2 = . Then substituting this value of p 2 
in (c) and integrating, we get 

(/) 4>t = (~ 3c, -f 0J)t - 2c 3 sin r - &0J sin 2r - T j T 0J sin 4t + c 4 . 

From the initial conditions <f> 2 must be zero when t = . Therefore c i = . 
Since it must also be periodic, c t = — \6\ . All of the constants of inte- 
gration are now determined except c, which will be determined by the perio- 
dicity condition on p t . 

The differential equations for p t and <f> 4 are the same as for p 2 and <f> 2 except 
in the right members. The process of integration is therefore the same. In 
the right members only even multiples of t occur except terms carrying the 
undetermined c 3 as a factor. In the equation corresponding to (i8d) there will 
be a term in cos t carrying c 3 as a factor. But the integration of this term will 
be non-periodic unless c 3 = . Put then c 3 = ; the integration proceeds 
just as before and the constants are determined in the same manner. This 
argument will be repeated in the coefficients of ft? and so on for all higher 
powers.* Therefore no odd multiples of t can occur in the solution. We have, 
therefore, 

p = 1 + [- i + i cos 2t + ^ cos 4t]0> 2 + • • •, 
(49) 

q = t + [- -/ T sin 2t - g|o sin 4t]0> 2 + • • • . 

* The differential equations arising at the successive steps are of the same type as the first 
two equations of (32), p. 555, in Professor Moulton's paper on periodic solutions of the problem 
of three bodies, these Transactions, vol. 7 (1906). Consequently his general formulas for 
the coefficients of the solutions, equations (42), loc. oit., apply here. 
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Since the series involve only even multiples of t the orbits are symmetrical 
with respect to both the r-axis and the g-axis. 

This completes the formal construction of the solutions of which the existence 
was proved in §§ 5, 6 and 7. 

Part II. 

ORBITS REENTRANT ONLY AFTER MANY REVOLUTIONS. 

§ 11. The differential equations. 

The orbits which we have previously considered have had the common prop- 
erty of involving only the period 27r . Since this period is independent of the 
oblateness of the spheroid the derivation of these orbits has been relatively 
simple. We shall proceed now to investigate a class of orbits which involves 
beside the period 27r another period 2irj\, where X is a function of the oblate- 
ness of the spheroid, the inclination of the orbit and the mean distance of the 
particle. We shall start from the solution which involves an arbitrary incli- 
nation. Into this solution four arbitrary constants were introduced, viz., incli- 
nation, mean distance, longitude of node and the epoch. Two more arbitraries 
are necessary for a complete solution, viz., constants corresponding to the eccen- 
tricity and to the longitude of perihelion from the node. In what follows we 
shall introduce the constant corresponding to the eccentricity, but we shall still 
project the particle from an apse at the node. 

We have found for the differential equations a certain solution (44) which we 
may write 



r = <f>(0, 1*;t), q = y}r(fi, /j,; r), C" ■■ ■■ c 



.2 _ „2 
0' 



which is symmetric with respect to the equatorial plane. That is to say, at r = 
the particle is in the equatorial plane and its motion is perpendicular to the 
radius vector. Its initial distance is <f>(0). Suppose now we change the initial 
distance slightly and also the initial velocity so that at t = 

r(P) = £(0) + a, ?(0) = 0, 

r'(0) = 0, ? '(0) = ^'(0) + 7 , 

and give an increment to the constant of areas so that c 3 = c\ + e . Can we 
determine these three constants, a , y and e , in such a manner that the series for 
r and q shall be periodic? These series can be expressed by 

»•= <K/3> m; t) + P, ?=|(^/*;t) + (t, c 2 = c 2 + e. 

If now we substitute these expressions in the differential equations (6) all the 
terms independent of p, a- and e will drop out, and there will remain the fol- 
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lowing differential equations for p and <r : 

(a) p"+{\ + [( - 40 2 + f/3 2 )- f/3 2 cos 2r] /4 2 + [(2O0f + 60 2 /3 2 + f/3*) 
+ (0 2 /3 2 - f/3*) cos 2t + f/3 1 cos 4t] /x 4 + • • • } p + { - 3/3/4 sin t 
+ [(_ 30 2 /3 + f/3 3 ) sin t - f/3 3 sin 8t] /t s + ■ ■ • } o- 
+ c {{ 1 + [( - 30 2 + 2/3 2 )- f/S 2 cos 2t] /4 2 + [(150*-i0 2 /3 2 + ff/3 1 ) 
+ (V-^yS 2 - -is/3 4 ) cos 2t + ^/3 4 cos 4r] m 4 + • • • } 
+ {-3+ [(120 2 -3/3 2 ) + 3/3 2 cos 2t]/4 2 + ... }/>+.••} 
= {3 + [( - 160 2 + 6/3 2 )- 12/3 2 cos 2t> 2 + [(9O0f + \9-0 2 /3 2 + |f/3 4 ) 

(50) + (i0 2 /3 2 - f f/S*) cos 2r + J e y./8* cos4r]/4 4 +. ■ • } p> + {-12/3/4 sin t 
+ [( - 900* £ + -Y-/3 3 ) sin t - -^-/S 3 sin 3t] p? + • ■ ■ } pa- 

+ (f + [(!*? - ¥ ^) + -V-/ 32 cos 2r] m 2 + • • • } ^ + 

(6) <r" + { 1 + [- f/3 2 + §/3 2 cos 2t]> 2 + [-40 2 /3 2 +70 2 /3 2 cos 2t]/4 j + • • ■ }a 

+ { - 3/3/4 sin t + [( - 30 2 /3 + f/3 3 ) sin t - f/3 3 sin 3t> 3 + •■■}/> 

= { - 6/8/* sin t + [( - 216> 2 /3 + \3-Z3 3 ) sin t - - 2 J-/3 3 sin 3t]/4 3 + . . -}/> 2 

+ {3+ [(30 2 - 3 4 3. /3 2 )+ ^_ /S 2 C0s2t -| m 2 + . . .} /3<7+{ | y8/iS ; nT+ .. . }<r 2 + 

In the first of these equations the coefficients of all terms containing odd 
powers of a- involve only odd powers of p and sines of odd multiples of t. All 
other coefficients involve only even powers of p and cosines of even multiples 
of t. In the second equation the coefficient of every odd power of a involves 
only even powers of p. and cosines of even multiples of t. All other coefficients 
involve only odd powers of p and sines of odd multiples of t. These properties 
play an important role throughout the entire discussion. 

§ 12. The equations of variation. 

Considering merely the linear terms of the differential equations (50) for p 
and a we have the equations of variation : 

(a) p" + { 1 + [ ( - 40 2 + f/3 2 ) - | /3 2 cos 2r] /4 2 

+ [(2O^+60 2 / 3 2 +-|/3 4 ) + (0 2 / 8 2 -| y 8 1 )cos2T+|^cos4T]^+ ...}/, 

(51) + {-3/3/t sin t+ [(-3<9 2 /3 -f- f/3 3 ) sin t-§/3 3 sin 3t]/x 3 + • ■ ■}*= 0. 
(6) <x" + { 1 + [- f/3 2 + f/3 2 cos 2t> 2 + ■ ■ • } a 

+ { - 3/3/4 sin t + [(- 3<9 2 /3 + f/3 3 ) sin t - f/3 3 sin 8t] p 3 + • • -}p = 0. 
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These equations admit a solution of the form * 

Pj = A.e^<j>.(T), 

a-. = A.e^ylr.( T ) (j = l, ■■■• 4) 

where \. is a root of the fundamental equation and where <£(t) and ~^(t) 
are periodic functions of t with the period 2tt. The four values of the \ (real 
or imaginary) are associated in pairs, equal in value but of opposite sign.f From 
the fact that t does not occur explicitly in the original equations (1) it is known 
a priori that one pair of the X ; has the value 0,J and consequently, if we choose 
the notation so that \ 3 = \ = 0, the two corresponding solutions will have the 
forms 

P»= A A*( T )> P«=^«l>4( T ) + 7, *a( T )]' 

We shall consider first the two solutions in which the X are not zero. Let 
us first substitute in (51) the forms 

P = e i ^<f>(T) (•=•-!), 

a = e' AT i/r (t). 
After dividing out the exponential there remains 

<f>" + 2i\<f>' + [1 - X 2 + a 2 fji- + a if i l + ■ ■ •] 4> + [a, /*+ a 3 ,u 3 + . . ■] V = 0, 
i|r " + 2i\ty' + [l-\ 2 +b 2 fi 2 +b i n i + -.]-f +[a lA i + a 3/ tt 3 + ...]<£ = 0, 

where 

a, = — 3/8 sin t , 

a 3 = (- 30 2 /3 + §/3 3 ) sin t - §/3 3 sin 3t, 

« 2 = (- 4 #> + I/ 32 ) - f/ 32 cos 2 ^ 

a 4 = (200* + 6<? 2 /3 2 + f/3 4 ) + (0 2 /3 2 - f/3 1 ) cos 2t + f/3 4 cos 4t, 

& 2 =-§/3 2 +§/3 2 cos2t, 

6 4 = a sum of cosines of even multiples of t. 

With respect to equations (52) it is known that <j> and yfr are periodic with the 

*Floquet, Annates Scientifiques de l'Eoole Normale Superieure, 2d series, vol. 12 
(1883), p. 47. 

t See Les MSthodes Nouvelles de la Mecanique Ctleste, Vol. 1, p. 193. 
tlbid., p. 187. 
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period 2tt, and that X vanishes with jx since the problem then reduces to the 
two-body problem in which the characteristic exponents are all zero. It can be 
shown that cf>, yjr and X are expansible as power series in /j, of the form 

j=0 j=0 i=l 

Substituting these expressions in (52) we find : 
Coefficient of ft . 

ft' + *. = 0, 



Therefore 

(53) 

Coefficient of ft. 

(54) 



$ = a^ cos t + a ( 2 0) sin t , 
yjr iy — 7^ cos t + 7 ( 2 0) sin t . 

#' + *i = - 2 *\ < / > u - Oi^o' 

-f r + -fi = - 2ix i'fo - «i^o- 



Since the periodicity conditions demand that the coefficients of cos t and sin t 
in the right members be zero, we must take X l — . We get then 

ft + 4> t = f/87?J - f/Sy^ooa 2t + f/S 7 f sin 2t, 
(55) 

^r" + yfr 1 = |/SaW - f /3af cos 2t + §/8af sin 2t. 

Integrating we have 

<^ = a™ cos t + a« sin t + |-0t«°> + i/3 7 f cos 2t - A/St"" sin 2t, 
(56) 

•f j = 7 < ,) cos t + 7 w sin t + f/84 ' + J/3a ( 2 °> cos 2t - -l/Saf sin'2T. 

Coefficient of fi 2 . 

,.„, </>' 2 ' + </> 2 = - 2 * ^o - «2 *« ~ «i -f i . 

^2 + ^2 = - 2i Wo -hfo- «1 4>1 ; 

or, expanded, 

ft + <j> 2 = B/Sy," - l/S^cos 2t + IV/) sin 2t + (40»af 2 °> + 2a 2 a< I 0) )sin t 

+ (4(9? a?" - 2tX<>) cos t + 3/3 2 a< u) sin 3t + 3/3 2 <»cos 3t, 
(58) V ' ' 2 2 ; 2 ' 

f 2 + t 2 = W + (V-^Tf + 2i\ 27 W)sin t + 2iX 2 7 2 °»cos t 
+ |/3a ( ;> sin 2t - s/S^'cos 2t. 



78 W. D. MACMILLAN: PERIODIC ORBITS [January 

In order to satisfy the periodicity condition we must have 

2i\ ap + 40J dp = , 2iX 2 yf + Y-/S 2 yf = , 

40 2 a ( 1 o >-2iX 2 a 2 °> = O, + 2i\ yf = . 

The equations of the second column are satisfied by taking 

7 co) =7 (o) =0 . 
Solving the other two we find 

X 2 =±26> 2 , «<"> - mf = . 

Equations (59) can also be satisfied by taking 

X 2 = dp = ap = yp = , -yf = arbitrary, 

but this would lead to the development of the solutions in which the characteris- 
tic exponent is zero, and these solutions will be discussed later. 

It was known at the outset that there were two values of X equal numerically 
but of opposite sign. We will choose the one with the positive sign. The solu- 
tion for the negative X can be derived from the solution for the positive X . The 
condition ap — idp = still leaves us with an arbitrary constant. Since the 
equations are linear this constant will enter the solution linearly and may there- 
fore be taken equal to unity. The arbitrary constant is restored in (79) after 
the solutions are completely developed. We will take then af> = 1 which makes 
dp = — i . Consequently 

(60) $ = cos t — i sin t , ifr = . 

Integrating (58) with these values, we get 

</> 3 = f/Sy" + af cos t + af sin t + J/SyW cos 2t - J/SyW sin 2t 

(61; + 3/3 2 cos 3t - 3/3 2 i sin 3t , 

yfr 2 = |/Sa 2 ') + yf cos t + yf sin t + |/3a« cos 2t - \$dp sin 2t . 

Coefficient of pK 

ti+^IW? + [- 2X 3 + ±6\(dp - idp)] cos r 

+ [2iX 3 +46> 2 (dp+idp)] sin T+|/Syf sin 2t— f/3yf cos 2t, 
(62) ^' + ^3= [f«?^-f^^-!i^ 3 ] + [-2a 27 (" + |#f 1 ')] cos t 

+ [2iX 2 7 ( 1 "+|/8 2 7 ( 2 1) ] sin t+ [_f «< 2 '/3+V-^ 2 /S+f J i/3 3 ] cos 2t 
+ [ f-af /8+ V 6> 2 /3- T 9 g/3 3 ] sin 2t- £ /8» ^ cos 3t- f /S*^ sin 3t . 
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From the periodicity conditions we must have 

- 2X 3 + 40 2 «> - W 2 ! >) = , f/3 2 ^ - 2i\^ = 0, 
2iX 3 + ±6\{i*f + «y>) = , 2i\ y v + |y8 2 7 m = . 

The last two equations can be satisfied only if 7^' = 7^ = . The first two can 
be satisfied only if X 3 = ( a ( p — ioQ>) = . The condition ( a ( ^ — ia ( ^) = again 
gives us an arbitrary constant. Then by (56) 4>i = c(cos t — i sin t), but this 
is the same as <f> u multiplied by c/m . That is, the solution is repeating itself one 
degree higher in n and this, of course, should be expected since the equations 
are linear, so that any solution multiplied by any power of fi must satisfy them. 
We are at liberty then to choose the arbitrary c = , which is the same as 
choosing a'^ = d^> = . Integrating (62) with these values, we find 

<£. ( = |/37f + af cos t + af sin t + ifiyf cos 2t — ifiryf sin 2t, 

(63) ^ 3 = [ 1^ _ |i0J /3 - f l iff ] + 7 f > cos t + 7 f sin t 

+ [J« ( 2 2 ^—V^i/8-t 7 6^ 3 ] cos 2t + C-K^-V-^^+i'e/ 33 ] sin 2t. 

It can be shown by induction at this point that <f> and X involve only even 
powers of fi, and that -ty is an odd series in fi. Furthermore 4> contains only 
odd multiples of t, and i/r only even multiples. Consequently 

7 f = 7® = a «» = af = 7 f = 7f = , 
and all 

Coefficient qfi* 4 . 

<f>' t '+<f> 4 = [_2X 4 +4^(af-kf)-160J-lO^/3 2 ] cos r+ [2zX 4 

(64) + 40f (iaf +af )+16i<9'] sin t+ [66> 2 /3 2 +3/3 2 af + T y3*]cos 3t 
+ [-6^ 2 /3 2 +3/3 2 of-ffz73 4 ]sin 3r-f i/3 4 cos5T + f l«/3 4 sin5T. 

Therefore, from the periodicity condition, we must have 

- 2X 4 + 40 2 ( of - iaf ) - 166> 4 - 1O0 2 /3 2 = , 

2 iX, + 46> 2 ( iaf + a» ) + lQi6\ = . 

Solving these equations, we find 

X 4 = - 80J - 1 0; /3 2 , af ) - iaf = |/3 2 . 

In this last equation we can choose dp = |/3 2 and a ( 2 21 = 0. This choice of 
o^ 2 ' and a ( 2 2) will make the coefficient of sin t in <f> 2 equal to zero, and since the same 
thing occurs for each <f>. it is evident that this method of choosing will simplify 



80 W. D. MACMILLAN : PERIODIC ORBITS [January 

the solution by making the coefficient of sin t equal to zero for all powers of fi. 
We have then on integrating 

<f> t = a t cos r + [- \6\{P - ^/S*] cos 3t + i [f 2 /3 2 + r %/3*] sin 3t 

(65) 

+ i|-g/3 4 cos 5t — T |^/3 4 sin 5t, 

*, = »"[-f^-!^] +*[~ : l 6 l ^-: t 7 6/3 3 ]cos2r+ [_i 6 J^/3-^/3']sin 2r, 

<f> 2 = f/3 2 cos t - f/3 2 cos 3t + §i/3 2 sin 3t, 

^ = - a i/3 - I i/3 cos 2t - i/3 sin 2t , 

O = cos t — i sin t , 

X= 26\^ + (- 80 4 - fW^ + • • -. 

The coefficient ot 4 of cos r in </> 4 is determined by the periodicity condition 
on (j> 6 . That the process of determining the values of the X. and the constants of 
integration arising at each step is general may be shown as follows. Let us 
suppose that we have computed everything up to and including </>. with the 
exception of the constants of integration in <f>.. We have then 

<f>. = dp cos t + dp sin t + known terms. 

The dp and dp enter i/r J+1 as follows: 

^J+i + tyj+i = 3/3 sin t <f> . + known terms, 

= §/3a« — |/8a ( /> cos 2t + §/3a«> sin 2t + known terms. 

Consequently in so far as it involves d/> and dp 

f. + i = |/3a«) + -J-/3a«> cos 2t - |/3^ sin 2t. 

Similarly, in so far as (f>. +2 depends upon constants as yet undetermined 

#+■+*,+»= - 2 ^;- 2i Wi + K 4 ^-!^ 2 ) + f£ 2 cos 2t J ^+ 3 ^ sin T t m 

= [- 2X j+2 + 46\ (d/> - M?) + A j+2 ] cos T 
(66) 

+ [2iX J+2 + 40 2 ( idp + dp ) + B j+t ] sin t 

+ 3/3 2 a ( /> cos 3t + 3/3 2 dp sin 3t, 

where A.,, and i?., „ are the known terms in the coefficients of cos t and sin t 
respectively. From the periodicity condition we must have 

- 2X J+2 + 40^ (dp - idp) + A J+2 = 0, 
(67) 

+ 2iX j+2 + \Q\ ( idp r + dp) + 5 = 0. 
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The solution of these equations is 

( 68 ) A -IB 

«(')_^) = _^_p«. 

As has already been pointed out we can choose dp = and we have then 

A^. - iB s 



(69) 



(./> zzi+1 j+ 2 



a\r = — 



80 2 



In order to show that X. +2 and dp are real it will be sufficient to show that 
A. +2 is real and B. +2 is a pure imaginary. This is readily shown by induction, 
for up toj = 4 inclusive we have 

0. = y^ ni K cos kt + * IC'** s ' n KT 

K K 

■xjr. = i £/, cos KT + Z •?« sin KT , 

where m K , »i K , f K and <7 K are all real. From the form of the differential equa- 
tions it follows at once that the same forms hold for j = 5, then j = 6 , and so 
on. That is A J+2 is real while B. +2 is a pure imaginary. 

It is further to be noticed that A +2 and B. +2 do not contain any terms in /3 
independent of 8\, and consequently the 8\ which appears in the denominator 
of dp will divide out. This is proved as follows. If 6\ be put equal to zero in 
the differential equations then equations (52) become the equations of variation 
of a circular orbit in the ordinary two-body problem, the plane of the circle 
being inclined to the plane of reference by an angle whose sine is /3/u. = s. 
The original differential equations (6) can then be written 

( 70) r " = ( 1 " -^~ ~ -- } - r -.-^v = S, <?"=-, - 2 -l tT ,= Q, 

where the constant c 2 is given the form (1 — s 2 ) ( 1 — e 2 ) . For these equations 
we have the solution 



(71) 



where 



(l_e 2 )yi-s 2 sin 2 (0-n) 
r _ a __ ___... ff _ j 

(l-e 2 )ssin(6>-fl) 
q = a 1 + e coMT-~0j~~ ' 



((?_5 ) = (T-T )+2e8m(T-T ) + 

Trans. Am. Math. Soc. 6 
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Let us now form the equations of variation by putting 

r = r + p, q = q + ff , e = e + e, 

where r g , q and e () are the values in (71). We find 

„ BR dR DR 

r dr r cq de 



„ BQ dR dQ 

or cq oe 



(72) 

Three solutions of these equations are given by 

(73) 



dr 


dr 


p ~° i aa ' 


*-c dq -» 
c ' oa ' 


2 5t„ ' 



5, 

c d Vo 
C »3e ' 



^e _ de _ 5e 

Ci sir e_C2 a To ' * s de " 



If e =J= these three solutions are distinct, but the case in which we are interested 
is when e = . In this case it is not difficult to see that the first two solutions 
coincide. Since the equations are linear, the system 

_ [ dr o <^o~l 
P ~ C *lda 5t J' 

n*\ „ , \ dq ° dq »l 

e - c *\da 5tJ' 

is also a solution, but as it vanishes for e = it carries e as a factor. We 
can divide out this factor and absorb it into the arbitrary c 4 . For e = this 
solution does not now vanish, and it is moreover distinct from the first solution. 
Thus we have three distinct solutions even when e = 0, but since BQ/de = 0, 
and dR/de Q carries e as a factor, equations (72) pass over to the equations of 
variation of a circle when e g = . For these equations we have three solutions 
which are periodic with the period 1-rr . The fourth solution is not periodic but 
involves a term of the form r times a periodic function. 

Let us return now to the solution which we have developed, (65), and consider- 
only the terms which belong to the two-body problem, viz., the terms which are 

*See Les Methodes Nouvelles de la Mecanique Celeste, vol. I, p. 163. 
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independent of 6\ . This solution may be separated into two solutions, one of 
which is real, the other a pure imaginary. The real solution is the third solu- 
tion of (73), and the purely imaginary is the second solution. Since both of 
these solutions are certainly periodic with the period 2?r , it follows that no term 
in /3 alone can occur in the A. +2 and B. +i of (66) for the presence of such 
terms would give rise to non-periodic terms in the two-body problem. Hence 
A. +2 and B. +2 carry &\ as a factor which may be divided out in equation (69). 
Further X. +2 , (68), carries &\ as a factor and therefore X vanishes with the 
oblateness of the spheroid. 

The solution (65) may now be written 

(75) 

where 

<£« = cos t + [ f/3 2 cos t - f/3 2 cos 3t] ^ 2 

+ [a 4 cos t + (- \6\P - fg^F) cos 3t + T | T /3 4 cos 5t]/^ + • • • , 

<p = - sin t + [|/3 3 sin 8t] fi 2 

(76) + [(f 0?/3 2 + JfajP) sin 3t - T | ¥ /3* sin 5t] /** + • • •, 

f o = [- 1/3 sin 2t] /* + [(- V-^x /8 - 1V/ 33 ) ^n 2t] /*« + • • ■ , 

!p>= [- f/3- 1/8cos2t]m 

+ [(- W " !^) + (- ii0 2 /3 - t V^)cos2t]^ + .... 
By putting 

e <\T __ cos ^ T _j_ j s j n Xt 

we can write 

pW = [<£<') cos x T — <p sin Xt] + i [<p cos Xt + <f> w sin Xt], 
(77) 

<7« =[yjrW cos Xt — i|r( 2) sin Xt ] + i [ ip cos Xt + f m sin Xt ] . 

We have thus one solution of the differential equations. A second solution 
can be derived from it by merely changing the sign of i. Thus 

f> = [<p cos Xt — <p sin Xt] — i \<j> c2> cos Xt + <f><» sin Xt] , 
(78) 

o< 2 > = [ ,p> cos Xt - ip sin Xt ] — i [ I// 2 ' cos Xt + ip> sin Xt ] . 

By adding and subtracting these two solutions we have finally 

P = A[ P ^ + P v]+ 5 [>»-/*»], 
(79) 
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A and B being two arbitrary constants. As above developed there is a certain 
arbitrariness in these solutions owing to the manner in which the constants of 
integration were determined. They may be reduced to a normal form by multi- 
plying each solution by the proper power series in y? with constant coefficients. 
By this process we can make for t = , 

^"(0) + /3 f2 )(0) = 1, o-W(O) - o- (2) (0) = fa. 

Since [/> (1) — /o (2) ] and [cr* + <7 (2) ] are sine series they vanish for t = . 

The third and fourth solutions of the equations of variation, (51), are given by * 

(80) 

«■ r dq « „ n d ( a V«) 

3 dr ' * da 

In forming the partial derivatives with respect to a in the fourth solution it 
should be remembered that 0- is an explicit function of a, by equations (4), 
and that t is also a function of a implicitly through n. The third solution also 
can be normalized by giving the arbitrary constant such a form that at t = 

ft = 0, * 3 = C0n. 

The fourth solution is non-periodic and has the form 

(81) /.-DJr^ + fc], «t« j d[t^»+*«], 

where (j> 4 and yfr i are periodic functions of t with the period 2ir. As in the 
previous solutions this may be normalized so that at t = 

Pi = DB*f, <r 4 =0. 

It is also easy to find (f> t and -ty- t by substituting (81) in the equations of varia- 
tion and solving for these variables, which must be periodic. 

Carrying out the above operations we find the following general solutions 

P=A{cos(1-X)t+[- : !-/3 2 cos(1-\)t+!/3 2 cos(1 + X)t-!/3 2 cos(3-\)t]ai 2 

+ [(-J« 4 + l^/3 2 +M/3 4 )cos(l-X)T+(ia 4 -||/3')cos(l + X)T 

+ (_ 30J/3 2 _ ^/3 4 ) cos (3 - \)t- t § 7 £*cob (3 + X)t 

+ 1 | 1? /^cos(5-\)t> 4 + ...} 

*See Les Methodes Nouvelles de la Mecanique Celeste, vol. 1, p. 163. 
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+ J B{isin(l-\)T+[(-i / S 2 -|^)sin(l-\)T- 1 5g/3 2 sin(l+\)T 

- &/8» sin (8 -X)t] /»»-..} 
+ C { [- J/3 2 sin 2t] m 2 + [ ( T V e\P - k& ) sin 2r + T V/3 4 sin 4t] M 4 + • . • } 
+ D\ {[IP + i£ 2 cos2r> 2 + [(f0;/3 2 - ^/S 4 ) 

+ ( _| 6\ /3 2 + iyS 4 ) cos 2t— s \^ cos 4t] m* • ■ ■ } +t{ [ f /3 4 sin 2t] /* 4 

+ [(!^?^ + |f^ 6 ) sin 2t - ^/3 6 sin 4t> 6 + . • . } }, 
o-=^{[|/Ssin\T- i/3sin(2 - \)t]/* + [§0 2 /3 sin\r 
--y-^ 2 /8sin(2-\)T]/ i 3 +...} 
+ R{ [|/8 cos \t + J/3 cos (2 — A.)t] /* +[- |0 2 /3 cos \t 

(83) + J<9 2 /3cos(2-\)t> 3 + ...} 
+ C{[/3cost>+[0]/* 3 +...} 

+ 2){{[i/8 sin T ]|t+ [(-|0f/3 + |/3 3 ) sin t> 3 + ■••} 

+ t{ [- f/3 3 cos r> 3 + [(- ff 0»/8" + !|/3 5 )cos T > 5 + •• . } }. 

§ 13. Non-homogeneous linear differential equations with periodic coefficients. 
If we have a set of homogeneous linear differential equations 

dx n 

(84) __.J = £^(; )a .., (i = l ,...,„), 

where the 0^. are periodic functions of t with the period 27T, we know from the 
writings of Floquet * that, if the roots of the fundamental equation are all dis- 
tinct, the solution has the form 

(85) x i = ±A j e^<f> iJ {t), 

where the A. are arbitrary constants, the a. are constants known as the charac- 
teristic exponents, and the $„ are periodic functions of t with the period 2ir. 

Such equations arise in dynamics whenever we study small variations from a 
known periodic solution. If we confine our attention entirely to the first powers 
of the variations the equations are linear and homogeneous, and are known as the 
" equations of variation." They have been studied extensively by Poincare in 
Les Methodes Nouvelles de la Mecanique Celeste. If the second and higher 

* Annates Scientifiquea de 1't.cole Normale Sup^rieure, 2d series, vol. 12 (1883), 
p. 47. 
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powers of the variation are considered equations of the same type arise but they 
are no longer homogeneous. It is necessary for our purpose to derive the 
character of the solutions of such equations when the non-homogeneous terms 
are periodic, even though the period be different from 2ir . We will suppose 
that the equations are the same as (84) with the addition of periodic terms. 

Case I. 

We will assume first that the period of the non-homogenous terms is 2tt. 
The equations are then 

dx n 

(86) ^/ = g^ + ^.(0 « = l,. ..,.), 

where r (t) and g { (t) are periodic with the period 2tt. The solutions of (86) 
may be written 

n 

(87) ^ = l4-^(0 + t(0 (» = l, ■••,»). 

If in the differential equations we change t into t + 2-w the equations remain 
unchanged, but the solutions become 

n 

(88) x. = X) A j ^^ 4>y ( * + 2tt ) + f. ( t + 2tt ) . 

Therefore equations (88) are also solutions of (86), and consequently 

d n 

jMt + 2tt) = E<^(* + 2tt) + g.(t); 

also 

d " 

Forming the difference of these equations we find 

(89) j t [+ i (t + 2K)-ir i (t)]=£d i .[+.(t + 27r)-+ J (t)] (. = 1,. ..,.). 

These equations are the same as (84) and their solutions have the same form as 
(85), that is 

n 

(90) ^(t + 2tt)- t,.(0 = T.Bj**WY 

3= 1 

The constants B. in this case are not arbitrary but depend upon the differential 
equations. They may or may not be zero. Equations (90) may be interpreted 
as showing that ^(£) is composed of two parts, a periodic part and a non- 



1910] ABOUT AN OBLATE SPHEROID 87 

periodic part. We may express it therefore in the form 

n 

(9i) ^(0 = <^(0 + Z^<M 0.4(0, 



J"=l 



where <o.(t) are periodic with the period 2tt, and the /;.(£) are functions which 
must be determined. Changing t into t + 2tt in equations (91) and forming 
the difference ijr.(t + 2tt) — yjf t (t), we find 



(92) ^{t + 2tt) - ^(0 = T,^^'Ut)[e^%(t + 2tt) -/<,(*)], 

= 22 -Sy e<v ' <£« ( f rom equations (90). 
Comparing the coefficients in these equalities we see that 

(93) e^f..(t+2^) -4(0 = 1, 

from which we can determine the character of the functions f i: (t). 
For this purpose let us define a new function X..(£) such that 

(94) V0 = «^(«) + nr^- 

Then by virtue of the relations (93) the \(t) are periodic with the period 2tt , 
for we have 

X„(* + 2tt) = eV-Jvf^t + 2tt) + f^-J^, 

= « a/ ./;(0 + 1 -_- 2a 7„, 

= \.(0- 

On solving (94) we find 

(95) /„(0 = «A(0 + ^:i. 

where the X.. are periodic with the period 2-rr . This expression for the f r (t) 
substituted in (91) gives 

(96) ir i (t) = ( o i {t) + ±B.4> i .X ij+ ±^^^ ii (t). 

The terms included under the last summation sign are merely terms of the com- 
plementary function. All the other terms are periodic with the period 2tt. 
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This form for the ^.(i) fails however if for any^' 

e 2 ^" = 1, 
that is, if 

a. = mod \/ — 1. 

In this event equations (93) for such values of j become 

(97) /„(« + 2w) -.4(0 = 1- 

We will define the corresponding value of X ; . by the relation 

(98) \,(t)=f ( .(t)-~. 

By virtue of (97) X.. is periodic with the period 2tt, for 

X ij (t+2^=f iJ (t + 2^)- t -±^, 

=4-(0- 2 ^ 

Solving (98) iorf..(t), we obtain 

(99) ,4(0=X,,.(0+ 2 ' 7r . 
Substituting this expression in (91), we find 

(100) i/<\(£) = periodic terms + ,, £ B . e^ $.. (t) , 

.) 

the summation in the last term to be extended over all j such that a. = 
mod V — 1 • 

We have then the general expression for the solution 

(101) x t = J2 A : ** 4> h - {t)+ periodic terms + ^- £ B j <& 4> { (t), 

j=i * j7r :i, 

where the summation in the last term is to be extended over all j such that 

a. = mod V ' — T. 

This result may be stated thus : 

Theorem I. If the d v {t) and g t (t) are periodic with the period 2tt and 
if the characteristic exponents a. are distinct and none of them congruent to 
zero mod y' — 1, then the particular solution is periodic icith the period 2tt. 
If the a. are distinct but some of them congruent to zero mod V — 1, then the 
particular solution may involve terms of the form t times the corresponding 
complementary function. 
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Case II. 

We will suppose that the g { (t) are periodic and that the period is different 
from 2tt . We will suppose further that g { ( t ) are expressible in the form 

(102) g t (t)= £ [o y cos (j+/3)t+b v sin (j+/3){] = e^^/^)+ e -^^/f («), 

where f ( P(t) and/\ 2) (£) are periodic with the period 2w and /3 is any real 
number not an integer. The differential equations may then be written 

dx n — 

(103) -^ = g ^.x, + «J= i «"/?)(0 + e->-Wf*>(t). 

Since the equations are linear we can consider the particular solutions depend- 
ing upon e' J ~' 1$t and e~" J ~' 1$t separately. The complete solution will be the sum 
of the two. Let us consider first 

dv n — 

(104) i = ^e iJ x J + e^'f^{t). 



dt , =] 



If we make the substitution 



x. = e^'^'y. 



equations (104) become after dividing out the exponential 

(105) d {l + i/=ly, = ±e.. yj +/<»(«)• 

Since 9.. and /X 1 ' are periodic with the period 2tt the discussion of the character 
of the y. reduces to Case I. The characteristic exponents of equations (105) 
are (a — V — 1/3), where the a are the characteristic exponents for the homo- 
geneous equations (84). Applying Theorem I we conclude that the y i are peri- 
odic with the period 2ir provided 



a. ^ V — 1/3 mod V — 1 (j = l, •••,»). 

If, however, for any value of j 

a. = i/— 1/3 mod V 7 — 1 

then, in general, non-periodic terms will arise just as in Case 1. 
The discussion for the second part of the differential equations 

dv n 

(106) -^=H0 ij x J + e- v -^ff(t) 

shows similarly that the particular integral is periodic if 



+ -l/^l/3 (j = l,...,n). 
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In equations arising in dynamics the characteristic exponents a. enter in pairs, 
equal but of opposite sign. Consequently for such equations this last condition 



«,+ — V— 1/3 modi/ — 1 u = l, •••,»), 

does not differ from 

a. 4=+i/ — 1/3 modi/— 1 0' = 1, • ■ •, n). 

We have then the following 

Theorem II. If the 6 {j are periodic with the period 2tt, and if the g ( (t) 
have the form 

9M = Z [°a cos (* + $)* + h i* sin (* + #)']' 

<mcZ j/* £Ae characteristic exponents a. are distinct and none of them congruent to 
± V — 1/3 mod V — 1 , then the particular solution has the form 

x.= ]£ [c im cos(m + 0)t + d. m sin (m + P)t~\. 

m 

If any of the a. are congruent to ± ]/ — 1/3 mod ]/ — 1, £/iew, £Ae particular 
solution will, in general, contain, in addition to periodic terms, terms of the 
form t times the corresponding complementary function. 

Let us suppose that one of the a. is congruent to ]/ — 1/3 mod V — 1 . 
Then, in general, the solution will contain non-periodic terms. If now the g { (t) 
contain also terms of the form J' ° (p + q/3)t, where p and a are integers, these 
terms will not, in general, give rise to non-periodic terms in the solution, for the 
a. will not, in general, be congruent to V — iq/3. Indeed, if /3 is not rational, 
for no integral value of q except unity will V — lq@ = «■• But if /3 (and 
therefore a. also) is rational, then certain values of q do exist for which 
a. = y'—iqft. Let us set 

a= .V — 1, /3= r , 

where i and j are integers relatively prime, and similarly for / and J. By 
hypothesis 



then 



■-. — -y. — mod 1 : 
J J 



— q -= =^ mod 1 



3 
unless q = 1 + U (I an integer). In this event a term of the form 
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Sli (p + <l@)t can also be written ^ (r + (S)t (r an integer) which involves only 
the first multiple of /3, and the solution in general involves non-periodic terms. 

Case III. 

When two of the characteristic exponents of the homogeneous set are equal 
the preceding arguments are not applicable. Let us suppose that the homo- 
geneous equations are 

dx " 
(107) -_i = g0 tf!By « = !,...,.), 

but that a n _ 1 = a n . The solution then has in general the form 

x i = n £c j e*S<i>.\t)+ [A+ Bt]e^<l> in (t)+ Be'^^it). 
Suppose now the given equations are 
dx. 



dt 



* = Z^*; + ^(0, 



i=i 



where the 6.. are the same as in (107) and the g i are periodic functions of t 
with the period 2tt. Let the particular solutions which depend on the g { {t) be 

Just as in Case I we find 

d n 

( 108 ) dt l+M + 27T) - ^(0] = X><, [fy(« + 27T) - *,(*)], 

and consequently 

(109) ^(* + 27r)-t i (0 = 2c;e^ i ,(0+ [i'+^je-'ilO+.B'e-'U*)' 

The terms included under the summation sign are obviously the same as those 
treated in Case I. Neglecting these we may write 

(110) ^ ( t ) = e-< </,,, (*)■/*.(*) + e-' • <£,,„_„ ( * )/ i(M _ 1} ( * ) + periodic terms, 

where f in (t) and J c "^ n _ l) (t) are functions whose forms are to be determined. 
Let us form the difference ^ ; (£ + 2tt) — "^.(O by means of (110) and compare 
the result with (109). From (110) we get 

(111) ^ + 2?r) ~ fi{t) = ^f^t + 2^) -/„(<)] ^'<U0 
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Comparing this with (109) we see that we must have 

^"fjt + 2,r) -f in (t) = A' + B't, 
(llz) 

Let us now define two new functions 

fJ . , , T A ' livB'e^ 1 , B'te^' 

(113) L V > J 

By virtue of the relations (112), X ia (i) and \ ( >i-i)(0 are periodic with the period 
27T , and consequently 

, , . ,, , . T 4' 27r J B'e 2 -" "I #'« 

.4(0 = e_ ""'\„(0 f e2 „-. _ i - /^_iy + e^-^1 ' 

(114) L V ; J 

/'., ,,(£) = e-'o'X., .-, + -= — r, 

where \ n (0 and \(„_,)(0 are periodic with the period 27r. Substituting these 
expressions in (110), we find 

f A' l-n-B'e^ B't "1 

^(0 = L^;._ 1) - (e ^_ i)2 + (e ^_ 1) Je"»^ i „(0 

(115) 

+ / e 2a;, I 3 1 - ) «' , "'0 ! („- 1 )(n + periodic terms. 

Comparing these terms with (109) we see that they are merely terms of the 
complementary function and that therefore 

x t = complementary function -f- periodic terms. 

These results hold provided a n ^ mod V — 1. In the event a n = mod \ — 1 
equations (112) become 

(116) f in {t + 2tt) -.4(0 = 4' + Z?'«, f i{n _^t + 2tt) -/<,_„(«) = B. 

It is necessary to give new definitions to the X functions. We will let 

(H7) MO =.4(0 - f, t 2 - A ~ 2 *~t, \ ( „-o(0 =.4-»(0 - £«, 

which give 

(H8) .4(0 = MO + f/ + - s/ 77 "^ /«.-«('■) = Mo(0 + 2^- 
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It is readily verified that \ a (t) and \ Un _ V) (t) are periodic by virtue of (116). 
Substituting these expressions in (110), we find 

(119) nt)=e a "' [[§,? + — i~' r *]«M*) +i 7r t ^n-v(t)\ + periodic terms. 



If we let 

A' - B'tv _ B 

we get 



^7T Z7T 



t,(0 = e °"'{ [C« + i^ 2 ] <k„(0 + ^</>,t„-i)(0} + periodic terms. 

This expression is not quite the same as t times the corresponding complementary 
function for in one term we have the coefficient JZ) instead of D. The theorem 
for this case then is 

Theorem III. If '0. and g.(t) are periodic with the period 2tt, and if two 
of the a are equal hut not congruent to zero mod V — 1» then the particular 
solution consists of terms periodic with the period 2tt plus a constant times 
the corresponding complementary function. But if two of the a. are equal and 
congruent to zero so that the corresponding part of the complementary function 
has the form 

then the particular solution consists of a periodic function plus a term of the 
form 

e^{(Ct + lD?)<t> in (t) + Dt4> i(n _^t)}. 

§ 14. Special theorems for the equations of variation. 

The foregoing theorems presuppose merely the conditions that the coefficients 
are periodic with the period 2ir. Further facts with regard to the solutions may 
be established when further facts are specified with regard to the coefficients of 
the differential equations. Our equations of variation (51) may be written 

where the notation with respect to the 0's has the following significance : even 
subscripts denote functions even in t, and odd subscripts denote functions odd 
in t ; one dash indicates that only odd multiples of t are involved, and two 
dashes indicate that only even multiples of t are involved. The solutions, 
equations (82) and (83), may be characterized in the same manner, and are 
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then 



(121) 



ft = A M T ) + B M T ) + Ga 3 (r) + Z)[« 4 (t) + to^t)], 
p 2 = A^(t) + B/3 2 (t)+C0 4 (t) + Z>[%(t) + t3 4 (t)], 
°-.= ^(t) + ^7 2 (t) + Cy 4 (r) + -Z>[7 3 (t) + t7 4 (t)], 

«T 2 =^ 2 (T) +Bl t (T) +CS 3 (T) +i>[a 4 (T) + T$(t)], 



where the notation is the same as for the 0's with the exception that in the 
first two solutions every integral multiple of t is increased by ± Xt, e. g. t 
cos (3 + X)t. On the3e terms the dashes refer only to the integral part of the 
coefficients of t. 

Suppose now we have the following non-homogeneous equations : 

(122) 

where ^(t) and/^r) are periodic with the period 2tt. Since the characteristic 
exponents are i/— IX, — l/ — IX, 0, 0, by Theorem III the solution has the 
form 

Pi = (Pi) + fi - (ft) + w i( T ) + « T «3 + &[!**«« + T «J' 

P 2 = (ft) + *, = (ft) + «,(T) + «^ 4 + Z>|> 2 £ 4 + t3,], 

(1J3) 

°"i=(°"i) + '?i = ( "i)+ W 3 ( T ) + «T74+ &[|T 2 y 4 +T7 3 ], 

*1 = K) + ^2 = K) + «4( T ) + OTS 3 + & [****» + ^,]> 

where (p ; ), (o" 4 ) indicate the complementary function and | ; and rj. are the par- 
ticular integrals, of which the (o. are the periodic parts. The a and b are con- 
stants which depend upon the differential equations. Since the %. and »? t . satisfy 
the differential equations we have 

^rfi(T) = f,(T), |r f, (t) = 7,f,(T) + ^(r) + sr(T), 
(124) 

^,(T) = ,,(T), ^ 2 (t) = ^(T) + ^(t) +/(T). 

Changing t into — t in these equations, we get 

^f,(-T)=— f,(-T), ^? 2 (-T)=-^ 2 | 1 (-T) + 5 3 ^ 1 (-T)-gr(-T), 

(125) J rf 

^^,(-T) — ^(-T), ^^(-T)=+F 3 f 1 (-T)-^, 1 (-T)-/(-T). 
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If now we make the additional hypothesis that g ( t ) is an even function of t 
and ,/'(t) is an odd function, equations (124) and (125) may be combined into 
the following set 

^[f.(T)-f 1 (-T)] = [f 2 (T) + f 2 (-T)], 



(126) 



^[^(T)+f,(-T)]=^[f 1 (T)-f 1 (-T)] + ^[, 1 (T)+, 1 (-T)], 

^-['?l( T ) + ^l(- T )]=['?2( T )-^(- T )]' 

^[^(T)-^(-T)] = « s [f 1 (T)-f 1 (-T)] + 5;[, 1 (T)+'» 1 (-T)]. 

These equations are the same as (120). Therefore 

^( T )-^(-^)=^« 2 ( T ) + 5 « 1 ( T )+C 5 3( T )+^[^( T )+ T «3( T )J' 
_ f,(T) + f,(-T) = 4/8 1 (T) + ^(T) + C^(T) + 2)[^(T)+T^(T)], 

(127) 

»/ 1 ( t ) + ').(- t ) = ^ i (t) + 57 2 (t)+C7 4 (t) + 5[7 3 (t)+t 74 (t)], 

Putting t = in these equations we get from the first and the fourth 

(128) = Aa 2 (0) + Da 4 (0), = A8 /t (0) + D8 t (0). 

Either A = Z> = 0, or the determinant ct 2 (0)S 4 (0) - S 2 (0)« 4 (0) = 0. But 
it is readily verified that the determinant is not zero. Therefore A = D = . 
If we suppose that £ 2 ( ) = ^ ( ) = (we shall be interested only in such 
cases), it follows from the second and third equations of (127) that 

= B%(0) +<7^(0), = B%(0) + C 7l (0), 

and hence B = C = . Consequently 

f ' (T) " ^~ T) = °' ^ (t) + ^ ( " T) = °' 

(lz») 

'?i( T ) + '?i(- -t) = 0, '? 2 (t)-t; 2 (-t)=0. 

We have then 

Theorem IV. If #(t) is an even function of t and f(t) is an odd func- 
tion ofr, and if £ 2 (0) = ? ? 1 (0) = ^> then ^(t) a?M^ V 2 (' T ) are even functions 
ofr, and £ 2 (t) a«(Z ^(t) are odd functions of t . 
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In the same way it can be shown that if </(t) is odd and f(r) is even, and if 
1^(0) = i? 2 (0) = 0, then f, and ij 2 are odd and £f 2 and r) 1 are even. 

Let us suppose now that g(r) is periodic and contains only even multiples of 
t, and that # /*(r) is periodic and contains only odd multiples of t. The general 
form of the solution will be the same as (123). £ , , f 2 , rj l , and r) i satisfy the 
differential equations 

^f.(T) = f,(r), ^f,(T) = \Ur) + OMr) + g(r), 

(130) 

^(r) = „,(t), J^ 2 (t) = ^(t) + ^.(t) +/(t). 

Let us denote | t .(T + 7r) by ^(t) and ij.(t + 7r) by »?](t). Then by chang- 
ing t into t -f 7r in (130) we have 

(131) 

J^;(t) = *?;(t), ^^(^--^.'(tj + ^c^-^t). 

From (130) and (131) it follows that 

(182) <* d 

The solutions of these equations, which have the form of (120), are 

£ - |; = ^a 2 (r) + 5a, (t) + C« s (r) + Z> [« 4 (t) + ra 3 (T)], 

| 2 - |; = ^(t) + #3,(t) + Cfi 4 (T) + D[B 3 (T) + t£(t)], 

(loo) _ 

r,, + t?; = ^ 7i (t) + 5 72 (t) + C7 4 (t) + Z> [y 3 (t) + t7 4 (t)], 

V,+ V' t = A\(t) + B\(r) + C\(t) + I) [? 4 (T) + tS 3 (t)]. 
Forming these expressions directly from (123), we get 

£i - I'i = »,(t) - o»,(t + w) - [_air + Iott 2 ] a 3 (r) - 6tt[t« 3 (t) + a/r)], 

& " fi = «/T) - «,(T + 7T) - [«7T + lW\$Sj) ~ 69t[t3/t) + £(t)] , 

(134) _ 

^i + Vi = <» s (t) + <b 3 (t + tt) — [ewr + 1&7T 2 ] 7 4 (t) — 6tt[t7 4 (t) + 7 3 ( T )] ' 

V2 + 7)' t = w 4 (t) + w 4 (t + tt) — [owr + J57t 2 ]8 3 (t) - 6tt[t5 3 (t) + S^(t)]. 
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Comparing (133) and (134) we see that 

.4 = ^=0, C= - [air + J&7T 2 ], Z»=-6tt, 

w i( t ) — <u i( T + i") = 0, « 3 (t) + <o 3 (t + 7r) = 0, 

<B 2( T ) — <B 2( T + 7r) = 0, « 4 (T) + (» 4 (T + 71") = 0. 

Therefore ^(t) and &> 2 (t) contain only even multiples of t while <» 3 (t) and 
&> 4 (t) contain only odd multiples, and by carrying this result into (123) we find 

£i = s i( T ) + OT s( T ) + b [i T2 « 3 ( T ) + T «4( T )L 

,„. ^ = S3 2 (T) + aT^(T)+6[^3 4 (T)+T^(T)], 

(135) 

^■= » 3 ( T ) + « T 7 4 (^) + 6 [j^(t) + t7 3 (t)], 

„, = s 4 (t) + «tS 3 (t) + 6 [|^ 3 (t) + tS 4 (t)] , 

and hence we have 

Theorem V. If g(r) contains only even multiples of "t andf{r) contains 
only odd multiples, then £, anc? £ 2 contain only even multiples of t and n l and 
r] 2 contain only odd multiples. 

If in addition to the above hypotheses we suppose that g (t) is an even func- 
tion of t and f(r) is an odd function, then f : and ?? 2 are even functions of t 
and £ 2 and n l are odd functions. Therefore b = 0. But if # (t) is an odd func- 
tion andy^r) is an even function, then ^ and n 2 are odd functions and £ 2 and 
t^j are even functions, so that in this event a = . 

In the same manner as above we prove 

Theorem VI. Ifg(r) contains only odd multiples of t andf(T~) contains 
only even multiples, then f , and £ 2 contain only odd multiples of t a?ic? 17, and 
v 2 contain only even multiples. Furthermore £ : , £ 2 , ^ anc? t; 2 are periodic 
with the period lir . 

If </(t) is of the form "YL.m, cos (J ± X)t and ,/*(t) has the form 
2j . ?i . sin (J ± X ) t then, since ± l/— 1 \ are the characteristic exponents of the 
homogeneous equations, the form of the solution is, by Theorem II, 

£i = Z)P { * 1,cos (£ ±x ) T + 2>f sin(&±X)r+ a (1 >Ta 1 (T) + a (2) Ta 2 (T), 
a- * 

& = l5?a»(i±X)T + Z^>in (*± X)t + a<'>T/3 2 (r) + «®t^(t), 

A A- 

»/, = 2>f cos (A ± X)t -f X''','' sin (&± X)t + Q! (,) t7 2 (t) + a^-ny^T.), 

A A 

t? 2 = £4; ) cos(&±X)t + £sf sin(£±X)T+ a^T^r) + « (2) ^2 ( T ) ^ 

A A 

but since g (t) is an even function and ./(t) is an odd function, £ 4 and ?? 2 are 

Trans. Am. Math. Soc. 7 
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even functions and £ 2 and 77, are odd functions. Therefore all the coefficients in 
(136) which have the upper index (2) are zero. But if #(t) were an odd func- 
tion of t and/'(T) an even function then all the coefficients in (136) which have 
the upper index (1) would be aero. Therefore 

Theorem VII. If g(r) is of the form 53 . m . cos (j ± X ) r and f(r) has 
the form £\ n - s ^ n {j ^ M T > where ± V — 1 X are the characteristic exponents 
of the homogeneous equations, then the particular solution has the form 

f 1 = J^p a cos ( a ± \ ) t + A-T'a v ( t ) , 

£2 = 2>» sin (° ± M t + a t@ 2 (t)' 
ii 

V l = J^P C sin (c ± \)t + At7 2 (t), 

C 

V 2 = Sp,/ C0S (^ ±X ) T + AtE^t). 

(i 

Also, 

Theorem VIII. If g(t) has the form "£. to. sin (j ± X)t and f("r) has 
tlieform ^.n. cos (.;± X)t, where ± V — IX are the characteristic exponents 
of the homogeneous equations, then the particular solution has the form 

£1 = I» n (« ± M T + Sra 2 (j), & = 2> 6 oob ( 6 ± X ) T + BtP^t), 

11, b 

v i = ^ l r c cos(c ±\)t + 2?t7,(t), ?7 2 = ^ r d sin (d ± X)t + Bt8 2 (t). 

c d 

It is understood that in the above two theorems a , 6 , c , <# , and j are integers. 

§ 15. Integration of the complete differential equations (50). 
It will be convenient hereafter to use the following notation for the solutions 
of the equations of variation : 

p = Aa 2 (r) + Ba^r) + C« 3 (t) + Z>[a 4 (r) + t<x s (t)], 

V " >!; ff =i7,(T) + ^ 2 (T)+ <7 74 (t) + Z>[7 3 (t)+t7 4 (t)]. 

The a.(r) and 7-(t) are characterized thus : 

<* 2 (t) involves only terms of the form cos [(2?i + 1)±X]t, 

" siu [2w±X] t, 

" sin [(2w+ 1)±X] t, 

" cos [2n± X] t, 

" sin [2w] t, 

" cos [2n + 1] t, 

li cos [2?i] t, 

" sin [ In + 1 ] t . 



7iO) 


44 


( 44 ii. 


«,(t) 


lb 


i a a 


7 2 0) 


.44 


4 u i; 


a s( T ) 


44 


(.4 44 


74( T ) 


it 


; 44 44 


a 4 (r) ' 


44 


4 44 44 


7s( T ) ' 


U 


44 44 *4 
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It will be convenient also to write the differential equations for p and <r as 

P " + e 2 p + 3 <x = e m e + e im e P + e mP 2 + e m p* + e w y +..., 

(138) 

a" + 4 <r+ 6 3 p = e i( y + 6 m pa + mO »+..., 

where all the 9's are periodic with the period 2tt, and 2 and 6 i contain only 
cosines of even multiples of r, and 8 3 contains only sines of odd multiples of t. 
On the right side of the first equation the coefficients of terms carrying odd 
powers of cr contain only sines of odd multiples of t. All the other coefficients 
contain only cosines of even multiples. In the second equation odd powers of 
<r have coefficients involving only cosines of even multiples. All other coeffi- 
cients contain only sines of odd multiples. 
The initial conditions are 

p = a, <r=0, |o' = 0, <r'=S. 

We will integrate equations (138) as power series in a, 8 and e, t entering into 
the coefficients. From Poincare's extension of Cauchy's theorem we know that 
these series are convergent for any arbitrarily chosen interval for t , = tS T, 
provided | a | , 1 8 1 and | e | are sufficiently small. The equations of variation 
involve the period 2ir/\. The solutions are not periodic unless \ is rational. 
Hence 9\ must be chosen in advance so that X is rational. We will suppose 
then that \ = JJK, where t/and _ZT are integers, relatively prime. Then three 
solutions of the equations of variation are periodic with the period 2Ktt. 

Since p and a are expansible in powers of a, S and e, we may write 

P = P m a + /Vo 8 + Pwi € + /V 2 + Puo a8 + P<rJ 2 + ftoi ae + PoJ e + Pw/ + • • '. 

°" = °-ioo a + °"o.o g + °"ooi e + "W** + °"no* 8 + °"o20 82 + °"ioi ae + o'onP 6 + <V 2 + ■ ■ ■ • 

Substituting these expressions in the differential equations (138) and equating 
the coefficients of similar powers, we have : 

Coefficient of a. . 

Piou + <Vioo + 6, 3°"ioo= °' 

°"'l'oO + ^4°".00 + 6 iPlW= - 

The solution of these equations is 

Pm = Mi?) + Ba A*) + Ca 3 (T) + D [a 4 (r) + ra 3 ( r)], 
"ioo = ^7x(t) + By t (r) + C%(t) + D [y 3 (r) + t 74 (t)]. 

In order to satisfy the initial conditions we must have, at r = , 
Pioo = 1 ' °"ioo = °> P'iw = ° ' °".'oo = ° • 
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From these conditions we find 

Aa 2 (0)+Da i (0)=l, 
Ba[(Q) + Ca;(0)=0, 
^7,(0)+ <7 74 (0) = 0, 
^7;(0) + D[7i(0)+7 4 (0)] = 0. 

The solution of these conditional equations is 



[January 



(139) 



(140) 



.6 = 0, 



C=0, 



A 7,;(Q) + 7 4 (Q) im 

■"- A -"MOO 5 



A = « 2 (0)[ 7 ;(0) + 7 4 (0)]-« 4 (0) 7 ;(0). 
Hence the solution is 

*,o. = ^SWxCt) + ^ ( ,?o [7 S (T) + ™ 4 (r)]. 
Coefficient of B . 
(142) 



(141) 



Po'lO + ^AllO + #3°"010 = , 
°"oio + ^4°"oiO •+" ^sftrtO = ° • 



These equations are the same as for the coefficient of a . From the initial con- 
ditions we must have at t = 



°"oio = ° > 



PolO = » 



/>010 = > 

The solutions are 

PolO = ^iW^) + ^S?0 K( T ) + ™ 3 ( T )L 
^0,0 = A$o%(T) + ^0 [7 S (T) + T7 4 (T)], 



= 1. 



(143) 
where 



A m = 

010 



« 4 (o; 

A 



A n - a A°J 

^010 ^ 



Coefficient of e. 
(144) 



Pooi + ^/V + <^°"ooi = ^oo 
°"'o'oi + 6, 4°"ooi + 3 Pmi = °' 



The right member, # 001 , is a periodic function of t with the period 27r. Further- 
more it involves only cosines of even multiples of t. Hence by Theorem V the 
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solution has the form 

(145) 

<r 00. = ^7 1 (T)+5 72 (T)+(77 4 (T) + i)[73(T) + T7 4 (T)]+ 75 (T) + aT7 4 (T), 

where a is a constant depending on d m , <* 5 (t) is a cosine series involving only 
even multiples of t , and 7 5 ( t ) involves only sines of odd multiples of t . 
From the initial conditions we must have at t = 

Pool = ° ' °"ooi = ° ' Pool = ° ' °"ooi = ° • 

Determining the constants of integration so as to satisfy these conditions we 
have the solution 



(146) 
where 



Pooi = ^o'oW T ) + ^SSi K( T ) + ™ 3 (r)] + a 6 (r), 
*oo. = ^i7,(t) + AZ [ 7s (t) + t7 4 (t)] + 7 6 (t), 



a e( T ) = a s( T ) - « a 4( T )' 7«(t) = 7 5 0) - «7 3 ( T )- 

It will be seen later that the value of <x 6 (t) for t = plays an important role, 
so that it is necessary for us to verify that it does not vanish. By hypothesis 
« 5 (t) is the periodic part of the particular solution for p in the differential equa- 
tions (144). Let us put in these equations 

Pooi = <£ + «™ 3 (t), o- 001 = 1 /r+ «T7 4 (r), 

where <\> and y]r are the periodic parts of the particular solution. We find 

</>" + 2 <f> + 3 ^ = - 2aa' i (T) + 1 + [(- 30? + |/3 2 ) - f-/3 2 cos 2t]/* 2 + . . ., 

f"+0 t ir + 6 3 <t> = -2a/3' 3 (T); 

or, 

<f>" + B t 4> + 6,f = 1 + [(- 30 2 + f/3 2 ) - 2(a + | )/3 2 cos 2t] /* 2 , 

^" + ^ + ^ = - 2a /3 sin t/x + • ■ -, 

where we have substituted 

a = a + a 2 p 2 + • • • . 
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The functions <f> and i|r can be expanded as series of the form 

<£ = 4>o + "M 2 + • • •, f = frU + fsV? + ■■■■ 

Then we get 

4>q + <£o = ! > so tnat <£ = 1 + c cos t, 

i/rj' + ^ = (3 - 2o )/8 sin t + f/3c sin 2t. 
Since yfr } is periodic we must have a = | , and then 

•«/r 1 = c, sin t — J/3c sin 2t , 
4*2 + 02 = — i c ^ 2 cos T + other terms. 
Since <f> 2 is periodic we must have c = . Therefore 

4> = a 5 (r) = 1 + ^P x {^), a = f + ^ 2 P 2 (^ 2 ), 

« t (r) = + [|/3 2 + J/3 2 cos 2t] m 2 + • • ■ • 
Consequently 
(149) a 6 (0) = 0,(0) - aa 4 (0) = 1 + /* 2 P 3 (^), 

an expression which does not vanish identically. 

Coefficient of a 2 . 

P200 + "2P2OO + "s ff 20U = ■"'200' 
°"200 + d * ff m + ^200 = ^200' 

The right members of these equations have the following expressions : 

^200 = ^ ( ,oo[^oo a 2 + 0no«27i + *o»tf] 

+ AZ A ZL 2 ^2oo « 2 (™ 3 + «*) + *„. { « 2 ( T 7 4 + 7 3 ) + 7,(t« 3 + a 4 ) } 
+ 2 ^7,(T7 4 + 7 3 )] 

+ ^[M™s + «,) 2 + <U™, + <M( T 7, + 7 3 ) + M^ + 7 3 ) 2 ] 
#200 = A &m a l + ^no a 2 7 1 + 2 020 7 2 ] 

+ ^MmC 2 ^^ + «J + 5 no I « 2 ( T 7 4 + 7 3 ) + Y,(™, + « 4 ) } 
+ 2S r w ,7 1 ('7 4 +7,)] 

+ ^[^( Ta 3 + « 4 ) 2 + C(™3 + «4)(^7 4 + 7 3 ) + #020 ( T 7 4 + 7 3 ) 2 ]- 
The initial conditions are 

P200 = ' ""200 = ' ^200 = ' ^OO = • 
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Since equations (150) are linear, and their left members have the same form as 
(51), the solutions will have the form 

Pm> = Aa 2 + Ba i + Ca z + D [to, + a 4 ] 

+ AtM?) + AZAftMr) + <>> 3 (r), 
(151) 

ff 200 = A 1x + B Vi + C % + D [ T 7 4 + 7,] 

Imposing the initial conditions we find 
B= C=0, 

t;(OK(O)-0,(O)[ 74 (o ) + 7 ^o)] , )2 

, *;(ok ( 0)-*,(0) [7,(0) + ?; (<>)] ,0, , (2) 

"r ^ -"-100 100 

. ^(° K( °)- < M Q ) [7 4 (0) + 7,;(0)] , (2)2 
1 (O ) 7 ;(O)-^(Q) a2 (O) 2 ^o)7 1 '(0)-^(Q)« 8 (0 ) , m . 

1J — 2i ^ioo ' ^ ^ioo^ 1 



(2) 

ioo"ioo 



a (O)7;(O)- ^(O)a,(O) , (2j2 
Substituting these values in (151) we have 

152 Pm = A[ ™ Xl(T) + A "° A "« X ^ + <^( T )' 

( * m = ^f oVl (T) + A%AZy 2 (r) + A^ Vz {t), 

where 



3' 



(T) = Mt) + wwg)z^[7i(0) + 7;(0)i^ (T) 



iv ■/ nv-y i A 



+ ti»w»-ww" w ,)^(, n , 
+ M° >*(°>-*;(°K(°> C7 , (T) + n , (T)] , 

and similar expressions for x 2 , y 2 and ai 3 , y z , which we shall find later [equation 
(172)] do not interest us. The characters of x x and y x are known with the 
exception of <f> l and \fr l which we will now investigate. The functions <£, and yfr l 
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are those portions of the solution of the differential equations (150) which depend 
upon the coefficients of A ( ll a . These coefficients are homogeneous of the second 
degree in « 2 (t) and 7,(t). In S. m and S m the expressions 

#,00' m and m contain only cosines of even multiples of t; 
5, m| , 6 im and 5 lm contain only sines of odd multiples of t; 
a 2 ( t ) has the form a 2 = ^ a k cos [( 2& + 1) ± X ] t ; 

7j ( t ) has the form y l == ]>^ ^ s ^ n [ 2& ± X ] t . 

k 

Consequently, so far as the coefficients of A^ are concerned, R wa and S. m have 
the form 

^200 = 1 £ < } «m 2fe + D «f cos [ 2& ± 2X ] t } ^ + • ■ • , 
£ 2 oo = (E^sin (2A + 1)t + £^sin [(2* + 1) ± 2X] t} A^ + ■■■. 

k k 

By Theorem II terms involving 2X give rise only to periodic terms in the solu- 
tion whose period is IKtt . ■ By Theorem V those parts of p m and a 2m depend- 
ing upon the terms in 2? 900 and S m which are independent of X have the form 

/ J =i ? l( T ) + CTa 3 ( T )< *=2>j(t) + CT7 4 (t) 

respectively, where p l and p 2 are periodic with the period 2tt. Consequently the 
characters of x x ( t ) and y l ( t ) are 

(153) » 1 (T) = P 1 (T) + Cl Ta s ( T ), 2/ 1 (t)=(? 2 (t) + c 1 t 74 (t), 

where ^(t) and (? 2 (t) are periodic with the period 2Ktt . 
Coefficient of a.8 . 

( 154 ) Piio + ^ftio + 3°"iio = -^uo' °"no + ^""no + ^aPno = #110' 
where 

*uo = 2 ^MV„ [« + ^no« 2 7, + m vU 

+ UToo<o + ^lo^o] [2^o„« 2 (T« 3 + O + tfnoM^a + O + SC^ + Ys)} 

+ 2^ %(T7 4 +7 3 )] 
+ 2A<» w Af> [^ (r« 3 + «J 2 + ^ 10 (t« 3 + « 4 )(t 74 + 7 3 ) + <U*74 + 7») 2 ], 
^ 110 = 2^^[# 200 a 2 2 + 5 1I0 a, 7l + 8 020 yU 

+ [^^oi ) o + ^^ ( o 1 1 ) o][2^o„« 2 (^3+«4) + ^ 1 „{7 1 (^ 3 +« 4 ) + a 2 (T7,+ 7 3 )} 

+ 29 020 7 1 (t7 4 +7 3 )] 
+ 2^ ( 1 2 o ) o^oi ) o[^ 2 oo(^3 + ^f + C(™ 3 + «J(T7 4 + 7 3 ) + C(^+ 7 3 ) 2 ]- 
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The functions JR U0 and S uo differ from i? 200 and /S^ only in the constants A. jk . 
The initial conditions impose the same conditional equations. Consequently the 
solutions differ only in the constants A.. k , so that we can express them at once 
without computation: 



(155) 



Pm = 24V,^>. (t) + [A$ AZ + A[l\A^ ]x 2 (r) + 2A™A™x 3 (t), 



where the s»;( T ) and 2/,( T ) are the same functions of r as before. 
Coefficient of 8 2 . 
By symmetry with the coefficient of a 2 it is seen that 



(156) 



/>020 = ^lfa(T) + ^V.^?0»,(T) + ^»k(T), 

"020 = ^^(t) + ^l ) .^? 1 ) .y 1 (T) + A^ y 3 (r). 



Coefficient of 'e 2 ; . 

Since the coefficients of the first powers of a and 8 were homogeneous in the 
A { of the first degree, the coefficients of a 2 , a8 and S 2 are homogeneous in the 
A.. k of the second degree. The coefficient of the first power of e is not homo- 
geneous in the A { . k , hence the second power is not homogeneous. But if the 
functions <* 6 (t) and y 6 (r) were zero the coefficient of the first power of e would 
be homogeneous. By symmetry therefore, we can at once write down the terms 
involving the A Vk to the second degree. To these must be added terms in the 
first degree and one term independent of the A„ k . 

The differential equations are 



(157) 



^002 + ^2^002 "I" "3°"002 — -"002' 

°"oo2 + e i a <m + e 3 Pm = S m 5 

-""ME = ^20o/'o01 • lloftlOl^OOl • 020^001 • 101 Pool ' 



"•" "o20°001 * 



C 3 .2 i S 

The terms involved in these expressions are shown in the following table : 

"002 



Term 


A™' 

^001 


4 m 4 m 

A OOl^OOl 


AW 2 
001 


A (1 » 
^001 


^001 


1 


Multiplied 
by 


Pool 


a 1 


2aj(™ s + a 4 ) 


( ™, + o 4 ) 2 


2a 2 o 6 


2o 6 (™s + «4) 


«: 


"200 


Pooi ^OOl 


a tfi 




(™3 + «4)(T 4 + r 3 ) 


7i" 6 
+ a i7s 


°6( T >'4 + 73) 


«6 7 6 


e uo 


ff Soi 


rl 


27i ( n t + y 3 ) 


(m + 7*) 2 


27,76 


2(T4 + 7s) 


7l 


"o20 


Pooi 








«2 


( ™ 3 + a 4 ) 


°6 


1O1 
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In order to obtain the ^ 002 it is necessary in the above table only to change the 
8... in the last column into 0... . 

The solutions of equations (157) may be expressed in the form 

,, CO ^=^l , ^( T ) + ^.^.^) + ^^) + ^.«>) + ^.*^)+ a!(| (T), 

(158) 

The coefficients of A^ 0J in the differential equations are homogeneous of the 
first degree in a 2 and 7 t , every term of which involves the first multiple of Xt . 
Hence the solutions for these terms by Theorem II involve non-periodic terms, 
and we can write 

x 4 (t) = P a ( T ) + c,to 1 (t) + c 3 t« 3 (t), 
(159) 

J^ 1 ") = *«( T ) + c i*7»( T ) + c 3 t7 4 (t), 

where (P 3 (t) and P 4 (t) are periodic with the period ZKir. As is seen from the 
table, x 6 (t) and y 6 (r) do not involve the X. They have therefore the form 

xJt) = «?.(t) + o.tol(t), 
(160) • 3k ' 

2/6( T ) = f? 6( T )+ C 4 T 74( T )« 

It will be verified in (172) that we do not need to know the character of x 6 and y b . 

Coefficient of ox . 

The differential equations are 



(161) 



P"\ oi + &2Pm + ^3°"ioi = Bm> 

°"i'oi + ^°"l01 + d iPm = ^101- 

^ioi = 02oo L 2 PiooPooi] + ^no[/'ioo "ooi + Pooi°"ioo] + 0o2o[ 2o "ioo°"ooi] + ^loiPioo' 

^JOl = ^200 [ 2 ^100 PoOll "+■ ^110 [Pl00°"o01 + /»001°'lOoJ + ^020 [ 2o "ioo°"ooi J • 

The following table shows the character of the terms entering into these 



expressions : 



R u 





■^ioo-^ooi 


aW aW +AW id) 
A 100 yl 001 ^ A \ oo A 001 


2^< 2 > -4< 2 > 
^100^*001 


^101) 


^100 


Multiplied 
by 


thw Pool 


<4 


200, ( TO, + Oi ) 


(!■«. + «4) 2 


Sc^Ole 


2o8 ( ra s + a 4 ) 


"200 


Pioo "ooi 

4" P 001 ff 100 


Oj7i 


7i ( r «3 + «4 ) 

+ <*, ( n't + r 3 ) 


(ra 3 +a 4 )(ry i + 7 s ) 


<*e/i4 a»7e 


a«(T74 + 7s) 


"no 


°100 "ooi 


y\ 


2/i ( 174 + n ) 


(V4 + 7 3 ) 2 


27i7e 


27 6 ( r 74 + 7 3 ) 


"020 


PlOO 








a. 


( T <*3 + «4 ) 


#101 



In order to obtain S m it is necessary only to change the iJk in the last column 
into 6 ijk . This table shows that R m and S m differ from i? 002 and S m2 only in 
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the constants A Vk . Since the initial conditions impose the same conditional 
equations as for the coefficient of e 2 we can at once write down the solution 

Pun = 2AZA% Xl (r) H^Z^Z + A&AZMt) 

+ 2^f» <) x 3 (r) + AZx^r) + AZx 5 (t), 
(162) 

"... = ^Z^M*)+ [^A% + Af m AZ ]y 2 (r) 

+ 2AZAZy 3 (r) + A%y,{T) + A^y^). 
Coefficient of he . 

This coefficient can be obtained by symmetry from the coefficient of ae by 
permutation of the first and second subscripts of the A.. k . Therefore 

Pon = ^A^x^t) + \_A<g t A% + A%AZ]x 2 {r) 

+ ZA&AV^iT) + AZx 4 (t) + AZx 5 (t), 
(163) 

"... = 2A% AV yi (T) + [A^AZ + Af l0 AZ ]y 2 (r) 

+ 2AZA%y a (r) + AZl/^r) + A^ y t {r). 

This concludes the computation of all terms up to the second order inclusive in 
a , 8 and e . In order to establish the existence of the periodic solutions it is not 
necessary to carry the computation further. 

§ 16. Existence of periodic orbits reentrant only after many revolutions. 

We have chosen the initial conditions so that at t = the particle is crossing 
the /o-axis orthogonally. It is obvious geometrically that if at any future time 
it again crosses the p-axis perpendicularly the orbit will be a closed one and the 
motion in it will be periodic. The conditions that the particle shall cross the 
/j-axis perpendicularly at t = T is that at this epoch 

p'=cr = 0. 

The equations of variation have the period 2_ff7r. Therefore we shall choose 
T = Kit . Since p is an even series in t, and a is an odd series, all the purely 
periodic terms in p and a are sines and consequently vanish at t = Kir . The 
terms which do not vanish must carry t as a factor. The conditions for 
periodicity give us two equations, namely, at t = Kir , 

P = ° = a m a + %io S + °W e + a m a2 + a no aB 

„ a4s + «02o 82 + a m ae + ^on^ + a oo2 e2 + ■ ' •' 

(164) 

°" = ° = ^uio 3 + 6 o,o S + Ki e + h m° 2 + b no aB 

+ &020 82 + b m ae + b ou & + b m e * + • • • ' 
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where a ijk and b i . k are the values derived from the series just computed. Their 
values are as follows : 

n — J(2) ,. n — AW 11 n — AW ti 

"100 — ■"■ 100 a ' "oio ~ ^010 " ' "ooi ~ ^001 a 1 

(165) 

b = A & v b = A (i) v b = A m v 

"100 — ^ioo " ' "oio — ^oio v ' "ooi — ^oni u ' 

n — A'-rfnr 4 AM AW v -I- AW 2 Z 
"200 — ^lOO^l ' "^100 "^100*2 ' ■ a - lOO^S' 

5 = ^0) 2 y 4. JU) ^(2) ^ 1 -d(2)2- 
"200 -"100^1 ^ 100 100^2 ' iOOi'3' 

a = 2A (l) A m x 4 T J_ (1) j! <2) -I- 4 <2) /I ' It + 2 J 4 (2) J_ (2) r 
^110 — •"^ 1 ioo^ 1 oio'^i ~ L-^ioo-^-oio 1 -^loo^oioJ *2 ~ ""-^100 ^010*3' 

6 = 2J. (I) ^l a) w 4- r A m A i2) 4- A (2) AM 1J+ 2J (2) J (2) v 
„ _ AM 2 ^ 4. AM AW 7 4- AW't;. 

020 -a - 010*1 ~ -"010 01 0*2 * -^010*3' 

A = ^(D 2 ^ 4. ^d) JU) £ 4. J.() 2 7/ 

°020 -^010^1 ~ ^010^010^2 ' ^"-OloyS' 

«,., = 2^ ^i*i + VAZA% + ^ ^]S 2 + 2^ ^ 2 > £ 3 + AM m x 4 + A%x„ 
K = *AM m AZy, + {AM m AZ + A%AZYy, + ZA^Af^ + A%y, + A%y„ 

«ou = 2AM AM mXi + [ A M AZ + Af l0 AM ]5j + 2^™ J®S S + ^.5, + ^ ( oTo^ 
6 M1 = 2^ AZ ~y x + [ ^V, AZ + J£> ^ ] t a + 2^ ^» y 3 + AZy 4 + A$ y* 

a — AM 2 z -l aw aw 7 4 /4 (2)2 -7 4 /1 (1) 7 4 /( l2) 7-4 7 

"002 ^"-001*1 I ^"OOl ^001*2 ' ^-001*3 ' ^001*4 ' ^001 » ^ *6' 

^002 = A^y, 4 ^o'Ji^S,^ + ^Si.% + ^001^ + A Zy 5 + y 6 , 
where 

u = Kir --, , v = a;7T7. , 
(166) 

Let us solve now the equation (164), p = , for e as a power series in a and 
8 . We obtain 

(167) e = e l0 a + e 01 S + e 20 « 2 + e n aS + e 02 S 2 + • • -, 

where the coefficients e r have the following values: 

_ a ioo 
"001 

_ _ a oio 



01 



e 20 



a ooi 

^200 "''OOI ^100 *101 ffi 001 j" ^002 100 



«\ 



001 



1910] ABOUT AN OBLATE SPHEROID 109 

_ _ a no a ooi ~ a oii a ioo g ooi ~ a ioi a oio a ooi + ^ g oo2 CT ioo a oio 

"ooi 

. ^W^OOl ^011 ^010 ^001 ' ^002*010 



a ooi 



Solving now equation (164), a = , for € in terms of a and S , we obtain 

(168) e = 6 10 « + e 01 S + e 20 a* + e u aB + i 02 S 2 + • • -, 

where the e {J have the same expressions in the b klm as the e have in the a klm . 
Subtracting (167) from (168), we have 

(169) 0= [i 10 -6 10 ]«+ [i 01 -6 01 ]S+ [i 20 -e 20 ]a 2 + [€,-€„] aB+ [e 02 -e 02 ]S 2 + . • .. 

We must now examine the coefficients of this series. 

a b A® u A (2) v 

r - -i _ ^ioo L \w _ _ loojf' _ fjioo _ fl 

L 6 10 e i0 J — a —J. ~ J(2) „ J(2) v — U ' 

, 1 _„, "ooi "ooi -^ooi"' "^001° 

r - _ -i _ %p ^olp _ ^oio M _ __oiq^ _ o 
L 6 01 C loJ ,, 7, — J (2) ,. J(2) „ 

"'ooi "ooi y± ooi"' -^-ooi y 

Both of the linear terms therefore vanish. The .computation of the second 
degree terms is somewhat more complicated. It simplifies matters somewhat 
if we observe that the e.. have the same expressions in v and the y. as the e ;j have 
in u and x r It is sufficient therefore to compute one and derive the other from 
it. From (165) we have 

_ a 2iio ^MO^lOl i ^002 a in 

20 ~ n a 1 o s 

"ooi ooi ooi 

Substituting in these fractions the values of the a.. lc from (165) we get 

a 200 _ _!_ [^(1)^(2)2^1 , Ad) A(1) A&'^l , ^(2) 2 J(2) 2 - 3 1 , 

a ~~ A(i? 100 ooi u ^ -^loo-^ioo^ooi u ^ -"-ioo^ooi u ' 
ooi -^-ooi *- ' 

a im a m _ 1 I n Am ACD Am a(%) '•'hi r_ A (2 > 2 A m A (2) —A m A &) A (i,2 l — 
_ _ _ — 4^l 100 ^l 100 ^ 001 ^i 001 +L ^100^001 ^ooi - a -ioo' a ioo- a -ooiJ u 

ooi ^.ooi L 

n 7 1 \ —2 A f2 » 2 A (2 > 2 ^S — A m A (2 > A (2} — — A (2)2 A {2) — 1 , 

V 1 ' 1 ) ^ Jl um J1 ooi u ^-ioo^ioo^ooi u ^-ioo^ooi u J' 

1 _ L [av? Aof^i , A&f jm a® - 2 + A^A^^ 
— , (2 f JL 100 -^ooi M + ^ ioo-^-ooi-^ooi M ^ -^-ioo^ooi M 



"002 _ 100 



, 4(2) 2 ^a) H* + vl® 2 ^ (2) ^^<2)2 X 6 

+ ^-ioo-^-ooi M ^ ^100 ^ooi u ^-100 M 



]■ 
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For the sum of these three expressions, there results 

i r x 

c 20 — /|(2) 3 I L 100 001 100 oou w 

(172) ° 01 

— r^4 (1) A {2) —A™ A m ~\A<® -4- vl <2)2 - 6 ], 

L 100 001 100 001 I 100 u ~ 100 u I' 

the coefficients of xju , x 3 /u and x 5 Ju vanishing identically. This is the reason 
why it was not necessary to compute x 2 , cc s , and x. in (152) and (158). 
Changing the x. into y { and u into v gives us — e 2u . Hence 



(173) 



I" e — e 1 = — - I [A m A (2) — A (2> A (,} 1 2 \- — — 1 

L 20 20j — . (2)3 j L 100 001 "^100 "^001 J w ^ I 



Hence 



- a^ \A m j <2) — a^ a w i r - - -i 4- ^ m2 r - — -1 1 . 

^MooL^ioo^ooi yl ioo yl ooiJ| u v I ^ ^ioo I u v \\ 
But [xju — Pi/v] and \_x 6 /u — yjv~\ both vanish since 

as is readily seen from (153) and (160). It is also seen on referring to (159) 
that \_x t /u — y 4 /v] does not vanish, but is equal to 

** _ h = c f 1 ^ _ Z?~| 

A (2) L4 (l) A {2) — A™ A m ~\ Vv v\ 

(. l 'V L 6 20 6 20j- ~ , (2 ,3 „ ^ • 

"^001 L J 

Without repeating the details of the computation we find in a similar way 

H 7t ,r- T ^oioL^ioo^ooi ^lOO^OOlJ ^^looL^oio ^ooi ^OlO^OOlJ -^4 i>4 I 

(17o) [e u -e u ] = -^ ---I, 

^001 L J 

■4© r Am Ad) _ A(2) Am -jr-, T. ~i 
itiR\ r- , n ^oioL^-oio^ooi ^oio^ooiJ ^ #4 
(176) [e 02 -e 02 ] = -^- 1 - - -I. 

^001 L -" 

Substituting these values in (169), we find that the second degree terms in a 
and 8 are factorable. Thus (169) becomes 

T ^ _ ¥±\ [A {2) a 4- A {2) 8 4- • • -1 l~( A (v > A i2) — A {t> A a) ) a 

/-IHH\ A(Zf\ u v L 100 ' ^oio u ' J LV^ 1 100 ^001 ^lCO^OOl/" 

7 + (A m A {2) — A ( -> A (1) )S-\ 1. 

T V-^-oio-^-ooi ^oio^ooi/ ' J 

There are therefore two real solutions for 8 as power series in a. Substituting 
these values of 8 in (168) we find the two corresponding values of e. In this 
manner we find : 



T tVsi 


{ solution. 




+ • 


■ ■ = a + • ■ 

a 5 — aa i 


«.(0) + 


+ ■ 


1 

. . = a + . . 

a. — act. 


• = ^o) a + 
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AW A (2> —A (i) A m 

£ 100 001 010_ 001 

— A m A& — A& A m 

010 001 010 001 

(178) 

A m A 1 ?' — A (i) A m 

100 010 1(10 010 

— A m A^ — A {2} AW 

J1 0I0 Wl 010 001 

Second solution. 
A c2) i' ( ) 

8--^- + ...-^ a+ ...- 7 ; ( o). a+ .., 

(179) ' 

ioo _ oio ^ i n /v i 

e = - ;<(27 « - -772)- °H = -a + ■ • ., 

001 001 

where a 6 and 7 6 are the quantities defined in (147), and 7,1(0) is the value of 
d^g/dr for t = . Thus one solution for e begins with the first power of a 
while the other certainly does not begin before the second, but in both solutions 
S begins with the first power of a. 

§17. Construction of the solutions having the period 2Ktt. 

We have just proved the existence of series for p, a and e proceeding in 
powers of the initial value of p (let us call this value e).* The series for p and 
a are periodic in t with the period 2Ktt, and since this condition holds for all 
values of e sufficiently small each coefficient is separately periodic. The series 
for p is even in t, and the series for a is odd in t. These series have the form 

P = Pi e + P 2 e 2 + p 3 e 3 + 

(180) a- = <r,e + (T 2 e 2 + <r 3 e 3 + • • •, 

e = e,e + e 2 e 2 + e 3 e 3 + •■ •. 

We may substitute these series in the differential equations and integrate the 
coefficients of each power of e step by step, and determine the constants in 
such a way that p and a shall be periodic, and satisfy the initial conditions 

(181) p(0)=e, <7(0) = 0, p'(0) = 0, <r'(0) = S, 

where B is a constant at present unknown but to be determined in the process. 
Substituting the series (180) in the differential equations (138), we find for the 
coefficient of the first power of e 

(182) Pi' + *,ft + Vi- = 0001*1, < + ^i + <Vi = 0. ■ 

* The reason for changing a to e is that this parameter corresponds to the eccentricity in the 
two-body problem. 
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Since by the condition of orthogonality p must be even in t and a odd in t, the 
solution complying with the condition is 

,1QON Pl = ^^^ + ^[^(^ + «,( T )] + «l[0^(T) + «.(T)], 

(183) 

<r, = -4 (1) 7i (t) + /> (,) [t 74 (t) + 7s (t)] + e^ary^r) + 7 ,(t)], 

where <* 5 (t) contains only cosines of even multiples of r and 7 5 (t) contains only 
sines of odd multiples of t. 

In order that this solution shall be periodic it is necessary and sufficient that 

Upon imposing this condition the solution (183) becomes [see eq. (146)] 

mo„ P l = A(i>a ^)+e l [a 5 (T)-aa i (T)] = AVa 2 (T)+e l a 6 (r), 
(184) 

<r, = ^ (1) 7,(t) + 6 t [7 5 (T) - 07 s (t)] = -4 (1) 7i (t) + e l74 ( T ). 

It remains to impose the initial condition that p { = 1 at t = . From this 
condition we get 

(185) ' 1 = A^a 2 + e~a 6 , 

where <5 2 and a 6 denote the values of these functions when r = . 
Coefficient of e 2 . 

Pi + e zPt + d i a 2 = ^001 6 2 + ^101 € lPl + <Wi + d mPi (T i + O 2o°1 = ^2' 

(186) _ 

< + 6 i °" 2 + 3 Pi = 6 2M Pl+ e noPl (r l + ^020 < = ^2 • 

Every term of R 2 and S 2 contains either A {1) , e l or e 2 as a factor. Arranged 
in this manner we have the expansions 

^2=^ (1,2 [« + 0n„«2 7 1 + 0O2otf] 

+ ^\ [ 2 ^200 «2« 6 + ^11.(7, «. + « 2 7 6 ) + 2^ om7i76 + 1O1 « 2 ] 

+ e i [^2oo a « + 0iio a 6 7. + ^o 20 7 6 2 + 0,o, «,] + e 2 [0 ool ], 

#2 = ^ <1)2 [0~2oo«2 + 0"no«27i + 0o 2 o7?] +-4 (1) e, [20 2 oo«2^+0no(7i« 6 +« 2 7 6 ) + 20 O2o7l76 ] 

+ e i [0 2 oo a c + 0,io « 6 7 6 + <?„2o7 6 2 ] • 

In order to understand the character of the solution of these equations we must 
examine the character of the various terms. The coefficient of A {1) in both 
H 2 and 8. 2 is homogeneous of the second degree in a 2 and 7l . Its expansion 
therefore involves terms carrying 2\t and terms independent of \. By The- 
orem II the solution for the terms in 2Xt is periodic. The terms independent 
of X are cosines of even multiples of r in i? 2 and sines of odd multiples of r in 
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S 2 . These terms have the same character as those in the coefficients of ef and 
e 2 and may be considered under the discussion of those terms. 

The coefficients of ^4 (,) e 1 in both R 2 and S 2 are homogeneous of the first 
degree in a 2 and 7l , all terms of which carry the first multiple of Xt. By The- 
orem VII the solution for p 2 will carry the term ra 1 ( r ), and for cr 2 the term 
T7 2 (t). Non-periodic terms of this character do not arise elsewhere in the 
solution, hence, in order to avoid them, we must take either A m = or e { = . 
If we choose A m = 0, then, by (185), e l is determined so that we must have 
€j = l/a 6 , thus agreeing with the first solution (178) of the existence proof. 
But if we choose e x = so that by (185) A m = l/3 2 we are in agreement with 
the second solution (179) of the existence proof. We will commence by devel- 
oping the first solution. This necessitates the choice A (l) = 0, e : = l/a 6 . 

First solution. 

Since A (V) = all terms in R 2 and S 2 which carry At or any multiple of it 
vanish. There remains 

E 2 = 6 ? [0«o a S + ^110 «.% + #020 To + 0i„i « 6 ] + «2^001 ' 

(187) S 2 =e\ [0 m al + 6 m a 6 y 6 + 0^] . 
We have also 

nasi o - a ^ T) a _%»( T ) e _.I 

(188) ^-~X~' ,_ 5, ' l ~5,' 

It is easy to characterize -ff 2 and # 2 . -ff 2 contains only cosines of even multi- 
ples of t, and S 2 contains only sines of odd multiples of t. Since p 2 is even 
in t and a 2 odd, the solution satisfying this condition is 

p 2 =^'a 2 (T) + Z»( 2 '[Ta3(T)+a < (T)] + [^(T) + « 2 Ta 3 (T)]+e 2 [a 5 (T)+aTa 3 (T)], 

<7 2 =^ 7l (T) + ZK 2 )[T 74 (T) + 7 3 (T)] + [? 2 (T) + a 2 T 74 (T)]+e 2 [ 7s (T) + «T 74 (T)]. 

In this solution the terms are grouped according to their origin. The first 
two terms are the complementary function. The third arises from the terms 
carrying ef as a factor. The fourth arises from the terms carrying e 2 as a 
factor. a. 2 is a constant depending upon the coefficients of ef in the differential 
equations. ot 5 ( T ) an( i 75( T ) are t ne same functions as in the coefficient of the 
first power of e. By Theorems IV and V, % 2 {t) and ? 2 (t) are periodic func- 
tions of t with the period 2tt and so constituted that V 2 ( T ) contains only cosines 
of even multiples of t, and ? 2 (t) contains only sines of odd multiples. 

In order that p 2 and cr 2 shall be periodic we must have 



IP> = _ a 2 -ae v 



Trans. Am. Math. Soc. 8 
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and this makes 

p 2 = A fV a 2 {r) + e 2 a 6 (r) + j? 2 (t) - « 2 « 4 (t), 

<t 2 = ^ 7i (t) + 6 276 (t) + r 2 (r) - a 2 7 3 (T). 

In order that we may satisfy the initial conditions we must have p 2 (0) = 0' 
and this determines e 2 so that 



A {i > is not yet determined, but it is obvious that it must be zero in order to satisfy 
the periodicity condition on the coefficient of e' in which the non-periodic parts 
arising from terms involving the first multiple of Xt carry A & as a factor. 
Hence the only way to avoid non-periodic terms of this character is to choose 
A {2> = . Anticipating this step then we have 

« 2 « 4 — n 2 , . , , . . v 

(191) , '_ . 

* 2 = ^V^7 6 (r) + ? 2 (r) - a 2 y 3{ r). 

Hence p 2 contains only cosines of even multiples of t, and & 2 contains only sines 
of odd multiples of t. 

It remains only to show that this process of integration can be carried on 
indefinitely. Assuming that up to and including p tl and <j\_j every p, and cr 
is periodic with the period 2tt and that the p. contain only cosines of even 
multiples of r and a, only sines of odd multiples of t, except that p t _ t contains 
the term A (i ~ l) a 2 ( t ) and <7._ 1 contains the term A (f ~ l) 7, ( t ) , it will be shown 
that the same conditions obtain for the next succeeding step. For p. and cr. we 
have from the differential equations (138) 

(192) + ^,,,(^7, + ',«,) + 2^^7,1 + *«, 

< + ^ + * S ft = ^ (i_1, [ 2 ^200^1 «2 + ^o(^7 1 + ^ I « 2 ) + 2^ 020 ^ 1 7 1 ] +%• 

From the properties of the differential equations it is readily seen that <J> { con- 
tains only known terms all of which are cosines of even multiples of t , and that 
M*. contains only known terms all of which are sines of odd multiples of t. The 
coefficients of A {i ~ l) are homogeneous of the first degree in a 2 and 7, , and con- 
sequently each term involves a first multiple of Xt . Their solution gives rise to 
non-periodic terms of the form Ta,(T) and t7 2 (t). They carry A ( *~^ as a 
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factor, and since terms of this type arise nowhere else in the solution we can 
make them disappear only by putting A {l ~ l) = . The solution for (192) then 
has the form 

n OQ ft=^ (i, *2(T) + ^ (i) |> 3 (T)+a 4 (T)] + foM + a/n^T)] + e^r) + ara^r)] , 
(193) 

«t 1 =^'J7 1 (t) + 2)W)[77 4 (t)+7 8 (t)] + [^{T) + aT % (T)] +6 < [7.(t) + ot7 4 (t)], 

where rj. ( t ) and f ( ( t ) are periodic with the period 2tt , and by Theorem V, rj. ( t ) 
contains only cosines of even multiples of t , and f ; ( t ) contains only sines of 
odd multiples of t. 

For p l and a. to be periodic it is necessary and sufficient that 

D {i) = — a. — <7e., 
which makes 

P;= ^ (i) <* 2 ( T ) + ^O) - «W T ) + e MT), 
(194) 

a- 4 =^w 7l (T)+Si(T)-a. 78 (T) + e.7 s (T). 

From the initial conditions we must have p { ( ) = . This condition deter- 
mines e. to be 

1 _ a. 

Thus the constants are uniquely determined. The /a ; and <r. have the properties 
assumed for those having lower subscripts, and the process of integration can be 
continued indefinitely. Every A u) is zero. Since no terms involving the Xt 
enter, the solution has the period Itt . But the orbits represented belong to 
the class of generating orbits from which we started. In other words, we set 
out with a generating orbit for which the initial distance was, let us say, r % , and 
we have found another generating orbit for which the initial distance is r -f e 
(e arbitrary). There is nothing surprising in this, for >* is a function of an arbi- 
trary constant /3 . Let us suppose we had started with a definite value of /3 , 
let us say y8 . This gives us a definite generating orbit with a definite initial 
distance r Q . Let us seek now the generating orbit for which the initial distance 
is r + e. If e is sufficiently small we can evidently give an increment e to /3 
which will increase r by the amount e . We have 

»•.=/(£.)> 

Expanding the right member of the last equation by Taylor's theorem we have 

Bf &f , 

- d/3 e + 2 d/3f + ' 
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and by inversion, since 



5/8 


(3=00 


e = - d je + ..., 
which we might write 


e = e l e + e 2 e 2 + • • • 


Then by substituting in the generating orbit 


/3 = /3„ + 


e x e + e 2 e 2 + 



and arranging as a power series in e , we would obtain the generating orbit which 
we sought. As these are the same conditions that are imposed when we seek 
new orbits through the equations of variation it is to be anticipated that one of 
the class of generating orbits satisfies the conditions. 

Second solution. 

We return now to equation (186) where the two methods of satisfying the 
periodicity condition presented themselves, and we shall continue with the second 
solution. We choose e, = , which we found determines A a> to be A m = l/a 2 . 
From (82) it is seen that a 2 = a 2 ( ) = 1 , and therefore A m = 1 . Hence in the 
second solution 
(195) ft = « 2 (T), ^ = 7!(t). 

Using these values of A m and e t , R 2 and S 2 of (186) become 

R 2 = { d m a l + ^iio^Ti + ^o tf] + e 2 e m , 
(196) 

^2 = [ #200 «2 + ^101 a 2 % + #020 ?1 ] • 

All of the terms in these expressions except e 2 9 m are of the second degree in 
a 2 and 7, . Therefore they involve terms carrying 2Xt and terms independent 
of Xt. m also is independent of Xt. In the solution the terms depending 
upon 2Xt are periodic by Theorem II. As for the terms independent of X, R 2 
contains only cosines of even multiples of t and S 2 contains only sines of odd 
multiples of t. These terms give rise to non-periodic terms in the solution 
which has the form 

p 2 = A^a 2 (T) + Z*»[t*,(t) + « 4 (t)] + </> 2 (X, t) 

+ K( T ) + «2 Ta 3 ( T )] + e 2 K( T ) + « Ta s( T )] ; 

(197) 

* 2 = A^ 7i (t) + Z> 2 >[ T 7 4 (t) + 7,(t)] + * S (X, t) 

+ [?,(t) + a 2 T7 4 (r)] + € 2 [ 7 .(t) + aT7 4 ( T )]- 
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In these equations <£ 2 (X, t) and ^(X, t) are the periodic terms involving X, 
r) 2 and f 2 are the periodic terms with the period 2tt ; a 2 the constant belonging 
to the non-periodic part, and the coefficients of e 2 are the solutions arising from 
the coefficient of e 2 in the differential equations. In order that this solution 
(197) shall be periodic it is necessary that 

i> 2) = — a 2 — ae 2 , 
and this reduces p 2 and <r 2 to 

(lyo) 

( r 2 = ^7 1 (T) + ^(X,T) + ? 2 (T)-a 273 (T) + e 276 (T). 

In order that we may satisfy the initial conditions we must have p 2 ( ) = . 
Hence, since a 2 ( ) = 1 , 

The constant e 2 is determined by the periodicity condition for the coefficient of e 3 . 
Coefficient of e 3 . 

(199) P :; + e 2?3 + e 3 * 3 = b 3 , <+ e^ + e tPa = s 3 , 

where 

^3 = Wi e » + m e 2Pl + 2 ^200/'l Pi + ^liofo/'l + "iP,] + 2 ^020 <7 l O "2 + 6 mP\ 

+ e 2mP\°l + e i2»P l °\ + ^030 °\> 
S S = ^200PlP2+~ d U0[°-2Pl + °-lfr] + 2 ^020 °-l°" 2 + ^SOotf 

+ ^lotf* 7 ! + 012O/>1°1 + ^030 »?■ 

In classifying the terms which belong to the expansion of R 3 and S 3 we bear in 
mind : 

First. The 6.. lc in i? 3 involve only cosines of even multiples of t, except those 
which are coefficients of odd powers of <r (i. e., j is odd). These involve only 
shies of odd multiples. The reverse is the case in the 8. Jk of S 3 . If j is even, 
0^ involves only sines of odd multiples of t. If j is odd, B. jk involve only cosines 
of even multiples of t. 

Second. p x and p 2 are cosine series. The terms independent of X involve only 
even multiples of t. 

<r, and a 2 are sine series. The terms independent of X involve only odd mul- 
tiples of T. 

It is seen then that among the terms in R 3 which are independent of X only 
cosines of even multiples of t enter ; and among the terms in S 3 which are inde- 
pendent of X only sines of odd multiples enter. In the process of integration 
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therefore two types of non-periodic terms arise, first, those coming from the 
terms which involve the first multiple of Xt, and second, those coming from the 
terms independent of X . It is important therefore to separate the various terms 
into three classes : 

(a) terms independent of X, 

(b) terms involving first multiple of Xt only, 

(c) terms involving multiples of Xt higher than the first. 

The solution for these last terms is periodic. We rewrite then the differential 
equations (199) as follows: 

p" a + 0,P s + 0,«r S = e 3^u + e 2.A( X ' T ) +/,(*, T )+/s( T ) + /A kX , T )> 

The coefficients of e 2 are, explicitly, 

/i( x . T ) = m a 2 + 2 ^ 2 oo a 2 a 6 + ^iinKX; + 7,«b) + 20 O2O 7,7 6 , 
5-i( x ^)= 23 2M a 2 a 6 + £ 110 (a 2 7 6 + 7 1 a 6 ) + 2£ ft20 7,7 6 . 

These terms are homogeneous in the first degree in a 2 and 7,, and conse- 
quently involve only terms which carry the first multiple of Xt . They are of 
importance since they carry the undetermined constant e 2 as a factor. The solu- 
tion for these terms has, by Theorem II, the form 

p = F l (\,T)+b l Ta l (r), <r = £,(X, t) + 6 1 t7 2 (t), 

where -^(X, t) and (^(X, t) are periodic and involve only terms carrying the 
first multiple of Xt. b x is a constant depending upon ,/"i(X, t) and ^(X, t). 
It is found, by a calculation not difficult, to have the value 

6, = 3 + p?9{p?). 

The functions / 2 (X, t) and <7 2 (X, t) have the same characters as ,/,(X, t) 
and ^[(X, t). They are considered separately since they are independent of e.,. 
Their solutions may be written 

p = F 2 (X, t) + & 1 to 1 (t), <t = £ 2 (X, t) + 6 2 t7 2 (t). 

The terms ,/* 3 (t) and ^ 3 (t) are independent of X. / 3 (t) contains only cosines 
of even multiples of t, while <7 3 (t) contains only sines of odd multiples of t. 
The solution for these terms has the form 

Finally, f^lcK, t) and g t (k\, t) involve only terms which carry multiples of 
Xt higher than the first. The solution for these terms is 

p = F i (kX,r), a=G 4 (k\,T). 
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The entire solution is then 

p 3 = A^a 2 (r) + D^[Ta 3 (r) + «,(t)] + e 3 [« 5 (r) + «t« 3 (t)] 

+ *2 W*, T) + ^(t)] + [2? (X, T) + ^(t)] 

+ [^(t) + 6 3 t«,(t)] + 2? (*X, t), 
(201) 

*,= ^ 3 ' 7i (t) + 2)W[77 4 (t) + 7s (t)] + e s [7 5 (t) + aT7 4 (r)] 

+ * 2 [#i(*. t) + ^(r)] + [# 2 (X, t) + 6 2 t 72 (t)] 

+ [G ! s (t) + 6 s 7 74 (t)] + © 4 (*X,t). 

All of the functions « : (t), 7 i (T), i^.(r), and G ( (t) are periodic. In order 
that p 3 and <r 3 shall be periodic it is necessary and sufficient that the coefficient 
of t« 3 ( t ) and T7 3 ( t ) , and the coefficient of Ta t ( t ) and T7 2 ( t ) be zero. That is, 

Z>< 3 >=-6 3 -«e 3 , 
and 

«--*■ 

by which condition the value of e 2 is determined. In order to satisfy the initial 
conditions we must have p 3 = at t = , and this determines A m , 

^) = 6 3 a 4 (0)-e 3 a 6 (0) + ^ 1 (0)- J P 2 (0)- J P 3 (0)- J P 4 (0). 

"1 

Thus all the constants are determined except e 3 , and the solution is 

p 3 =^' 3 )a 2 (T)-6 3 a 4 (T)+e 3 « 6 (T)-^ 1 (X, t) + F 2 (\, t) + j P 3 (t)+ j P 4 (^X, t), 
(202) ' 

S=^ (3) 7 1 (t)-& 3 7 3 (t) + 6 3 7 6 (t)-^ ^(X , t)+ G 2 (\, t) + G s (t) + ##X, t) . 

The constant e 3 will be determined in satisfying the periodicity condition for 
the coefficient of e 4 . This process of integration can be continued indefinitely. 
p 3 and <r 3 have the same properties that have been stated for p 2 and a 2 . It is 
evident from the properties of the differential equations that these properties 
persist for p i and <r 4 , and so on indefinitely. The coefficient for e i _ l in so far 
as it carries the first multiples of Xt is always the same as for e 2 . Therefore 
the constant e i _ l can always be determined so as to avoid non-periodic terms of 
the type ra, ( t ) and T7 2 ( t ) . The constant D (i) of integration can always be 
determined so as to destroy non-periodic terms of the type ra 3 ( t ) and T7 4 ( t ) . 
The constant A^ can always be determined so as to satisfy the initial condi- 
tions. The analysis of the type of terms entering is the same as for the sub- 
script 3. 
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We have therefore a periodic solution with the period 2AV which is different 

from the class of generating orbits from which we set out, for the particle makes 

many revolutions before its orbit reenters. After integrating the equation 

dv/dr = c/r 2 the solution will contain five arbitrary constants (subject to the 

condition that X must be rational), corresponding to mean distance, eccentricity, 

inclination, node, and epoch. 

University of Chicago, 
April 20, 1909. 



